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PREFACE 

THIS book contains the substance of a course of lectures 
delivered to engineers and physicists at the National 
Standards Laboratory, Sydney, in 1944, when the author 
was seconded to the Radiophysics Laboratory of the 
Council for Scientific and Industrial Research. A mimeo- 
graphed edition was produced by the C.S.I.R. in 1945 
for circulation in Australia. 

There is an increasing demand for a knowledge of 
operational methods by engineers and physicists who 
encounter transient problems, but whose mathematical 
equipment is restricted to that of an ordinary engineering 
course. Of the many books available, the older ones use 
the rather suspect and certainly obscure " operational 
methods ", while the newer ones require a knowledge of 
the theory of functions of a complex variable, and some 
comparatively advanced mathematics. Here I have at- 
tempted to go as far as possible using no mathematics 
beyond ordinary calculus, but at the same time catering for 
the reader who wishes to acquire the manipulative skill 
necessary to solve his own problems : to this end, the book 
contains as little theory as possible ; it is, in fact, largely a 
collection of worked examples illustrating the methods of 
solution of the various types of problem commonly arising 
in circuit theory ; a handful of problems, with answers, 
has been added to enable the reader to test his skill. 

For problems on ordinary linear differential equations, 
and thus on electric circuits with lumped constants, the 
use of the theory of functions of a complex variable is 
quite unnecessary ; but it is needed for a complete study 
of partial differential equations, and in the higher parts 
of the subject. The engineer who wishes to go deeply 
into these matters must study this theory, and since quite 
sophisticated points in it do arise, he must study it fully 
and carefully, and so before he starts he is entitled to 
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know what types of problem he will be able to solve, and 
how useful the solutions will be. I have therefore ^ in 
Chapter IV discussed the Equations of the transmission 
line from the point of view of the earlier chapters, using 
calculus only ; this treatment covers all the common, 
problems, and will enable the reader to follow published 
work on the subject. It is to be understood that solutions 
obtained in this way are correct, and that the algebra is 
the same as it would be if the complex variable were used 
the latter merely changes the point of view, and allows ^ 
complete mathematical discussion of the validity of the 
solutions. 

The reasons for the notation used here for the Laplace 
transform roo v 

e-**y(t)dt . . (*) 

Jo 

of y(t) may be given briefly. Any notation should satisfy 
three conditions : (i) it must be short, and easy to write ; 
(ii) it must be applicable to any symbols ; (iii) it should 
not be used in any other connection. The simplest 
notation of all'would be to write a capital letter Y for the 
Laplace transform of the corresponding small letter y ; 
but this may be most inconvenient since it is not applicable 
to many symbols, in particular to capital and many Greek 
letters. A " star ", y*, is not much used in other con- 
nections but is rather clumsy to write. The " bar 
notation, y, used here, was originally used by Le"vy and 
reintrodueed by Carslaw ; it is easy to write, and can 
be applied to any symbols, but it has the disadvantage of 
being employed in other senses, these, however, rarely 
occur in those fields where the Laplace transformation 
is most useful, and ease of writing is a most important: 
desideratum. If preferred it can be replaced by y used 
by Millman. 

A much more difficult question is that of the relative 
merits of the form (i) and of the "_ -multiplied " form 

fOO 
e- pt v(t}dt 
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Both notations have their own minor advantages, but the 
essential point is that form (2) gives an exact correspondence 
with the large body of work done by the older operational 
methods, while (i) corresponds with the newer books, 
with most of the theoretical work, and with corresponding 
work on the Fourier transform. 

Finally, it is a pleasure to acknowledge the assistance 
given to me by many officers of the C.S.I.R., in particular 
that of the Secretary in making available the drawings 
for the text-figures. 
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CHAPTER I 

FUNDAMENTAL THEORY 

i. The Laplace transform 

In this book we shall primarily be concerned with 
problems on the behaviour of 'dynamical or electrical 
systems which are started in a given way at some instant, 
which we take to be the zero of time, t = o. Mathemati- 
cally, these are called " initial value problems ", and may be 
expressed as the solution of a system of differential equa- 
tions with given initial values at the instant t o. Negative 
values of the time do not arise in either the specification of 
the problems or in their solution. The Laplace trans- 
formation is a mathematical device which is useful for 
solving such problems : if a function is specified for all 
positive values of the time, we can (in the simple cases 
needed here) write down a related quantity called its 
Laplace transform ; conversely, if the Laplace transform 
is known, the value of the function can be found from it. 
The importance of the method arises from the fact that, 
for the types of problem we are interested in, it is easy to 
write down the Laplace transform of the solution, and 
from this the solution itself is then easily found. 

We suppose y(f) to be a known function of t for values 
of f> o. Then the Laplace transform y(p) of y(t] is 
defined as 

f*rf\ 

. . (i.i) 



poo 

y(p) = 

Jo 



where p is a number * sufficiently large to make the integral 
(i.i) convergent. For .example, if y(t) i any positive 

* Of course p may be a complex number whose real part is 
sufficiently large to make (i.i) convergent, but, in the early parts 
of the theory, it is a little more definite to think of p as a real 
positive number, and nothing is lost by this restriction. The 
Laplace transform is one of a number of " integral transforms " 
formed as in (i.i) by multiplying a function of t by a chosen 
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p is permissible, but if y(t] = e 2t we must have p > 2* 
Apart from this p is unrestricted, and y(p] is a function 
of just as y(f) is of t. We shall occasionally write them 
in full in this way to emphasize this dependence, but 
usually shall write y for the Laplace transform of y, it 
being understood that these. are functions of p and t 
respectively. 

We need first a collection of the Laplace transforms 
of a few common functions. For example, if y = eS 
where a may be real or complex, 



- - , (1.2) 
p - a 

where p must be greater than the real part of a. 
Again, if a is real and 

y = sin at = .(e iat ~e- iai ), 
using (1.2) we have 






2i\p ia p + iaj p 2 + a* v JX 

Proceeding in this way, or quoting the results from 
known definite integrals, we can build up the Table I 
of Laplace transforms which contains all the -transforms 
needed in the first three chapters. This Table, and 
Theorems I to III below, provide the whole of the theory 
needed in Chapters I and II. 

In (1.6), as remarked above, p has to be greater than 
the real part of a, while in (1.7) to (1.12) a is real, and in 
(1.9) and (i.io) p must be greater than* \a\. All the 
results of Table I are elementary integrals except (i.il) 
and (1.12), which have been included because they arise 
in cases of resonance. The simplest way of deriving 

function of t and p, and integrating with respect to t. The 
Fourier transform (26.1), is another example. Each integral 
transform has characteristic properties, but those of the Laplace 
transform which are proved below make it the moat suitable for 
the solution of initial value problems. 

* If a is real, \a\ a, if a > o ; |a| = a, if a < o. 
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results of this type is as folio 

states that 

poo 

1 e~ vt cos at 
Jo 

TAB 
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ws : written out in full (1.8) 


LE I 


y 


y(p} 


i 

t n-l 


J <'*> 


e at 

sin at 
cos at 
sinh at 
cosh at 

sin at 
za 

5 (sin at at cos at} 


- 1 . . (i fc\ 


p-a ' ' (I>6) 
CT -,\ 


P / T o\ 


a CT n\ 




p- a- v 

P Cl IT^ 


(p* _|_ a 'i)a t 1 - 11 / 
, T / T T -\ 


(pi _(_ ^aja U-is; 



If we differentiate * both sides of t this with respect to 
as if a were a variable, we get 



* This technique of differentiating under the integral sign 
with respect to a parameter is a powerful weapon of pure mathe- 
matics ; naturally its use needs justification, and this can usually 
be supplied when, as here, it is employed in deducing new Laplace 
transforms from known ones. 
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which is equivalent to (i.n). Treating (1.7) in the same 
way gives (1.12). 

The collection of transforms in Table I can be greatly 
extended by the use of the following Theorem. 

THEOREM I. If y(p) is the Laplace transform of 
and a is any number, real or complex, then y(p + a) is 
Laplace transform of e~ at y(t). 

This follows immediately, since the Laplace transform 
of e~ at y(t) is 



f V* . e~ at y(t}dt = \ er<*+*)t y (t)dt = y(p + a). 
Jo Jo 

As an example of the use of this Theorem, we see from. 
(1.5) that the transform of 



(- !)!"(# + fl)' 

Similarly from Theorem I and (1.8) it follows that the 
transform of 

e~ vt cos at is , . , N2 s, 
(p + b) z + a 2 ' 

and from Theorem I and (1.7) that the transform of 

. . , . a 

e~ u sin at is ... 5. 

O + *) + 2 

2. To j^c? the function y(t) which has a given Laplace 
transform. 

It will be seen in 4 that the Laplace transformation, 
method of solving an ordinary linear differential equation 
with constant coefficients and given initial conditions 
consists of writing down an algebraic equation for the 
transform y(p] of the solution y(t}. This transform will 
usually be found to have the form 

M (ail 

g(t). ' ( 2 -*> 

where f(p] and g(p) are polynomials in p which have no 
common factor, the degree off(p) being lower than that of 
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Most of the labour * of the solution consists of finding 
the y(t) corresponding to a y(p) given in this form. This 
is done by looking up y(p) in the Table I of transforms, 
possibly combined with Theorem I. It is shown at the 
end of this secton that the y(t) found in this way is the 
only function which has y(p] for transform. 

For example, if 

y = ; , it follows from (1.6) that y = e~ 2t . 

p + 2 v ' J - 

Similarly, from (1.7) and (1.8), if 

2* -f- 6 
y = p* + A> y = z cos 2t + 3 sm 2t. 

If the denominator of y is a general quadratic in p, 
we complete the square, and use Theorem I and the Table. 
For example, if 

# P + 7 P + 7 j> + i 6 



it follows from Theorem I and the Table that 
y = e~* cos it -f- 3e~* sin 2Jf. 

Notice that both numerator and denominator of y have to 
be expressed in terms of p -f- a (in the case above, p -\- i) 
before Theorem I is used. 

If the denominator of y is of higher degree than the 
second, we have first to find its factors,f then we express 
y in partial fractions by the usual algebraic methods, 
and finally write down the function of t corresponding 
to each fraction using Table I. Alternatively, Theorem II 
below may be used when it is applicable. 

* The student is particularly advised to make himself skilful 
in the processes of this section. Labour can be saved by the 
use of a more extensive table of transforms than that given here, 
but transforms not included in the tables at present available 
frequently arise, and in any case a solid background of manipulative 
skill is extremely useful. 

t For the method of doing this in problems where they are 
not obvious see 12. 
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Assume * 
Assume 



Then P+-IS* A(p z + 4^ + 8) + p(Bp + C) . (2.2) 
Equating coefficients of p^, p, and the constant term, 
respectively, we have 

A + B = o, 

8A = i! 

Solving these for A, B, and C, we find 
P-4 



- (2-3) 
~^j 

Then from Theorem I and the Table, we find 

y = e~ 2i {cos 22 3 sin 2t}. 

The process of equating coefficients used above is 
always available, but the results can often be obtained more 
quickly otherwise. For example, if we put p o in the 
identity (2.2) we find A = 1/8. Then, putting A = 1/8 
in (2.2), this becomes Bp + C = (4 p), and we get the 
result (3.3) obtained before. 



y = 



(p + + 4) . 



Then i m 

Equating coefficients we find 



* It is shown in text-books on algebra that if we assume a 
fraction with a constant numerator corresponding to each linear 
factor of the denominator, and a fraction with a linear numerator 
corresponding to each quadratic factor (or squared linear factor) 
of the denominator, and so on, we shall always have just the right 
number of equations to determine the unknown constants in 
these fractions. 
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Then, using the Table and Theorem I, we get 

y ~ %-* + %te~* $ cos 2t -^ sin 2i. 

There is a simple general Theorem applicable to many 
important cases : 

THEOREM II. If y(p} = f(p)/g(p), where f(p) and g(p) 
are polynomials in p, the degree off(p) being less than that of 
%(p), and if 

(P) ($ ^(P a z) - - (P a n \ - (24) 
where a^ a z , . . ., a n are constants, which may be real or 
complex but must all be different, then 



" ^- 

= Y L_ P 

^I^- fl r' L 



(P ~ 



. . . . (2.7) 
the function y(t) whose Laplace transform is y(p) is 

S-5?f\ eV (2-8) 

\ T ^ (2-9) 



(2.8) is the form taken by Heaviside's Expansion 
Theorem in the present method.* Since (2.8) . . . (2.10) 
follow directly from (2.5) . . . (2.7), respectively, by 

* In Heaviside's notation, and in the "^-multiplied " form of 
the Laplace transformation, in which (in order to have an exact 
correspondence with the Heaviside notation) the Laplace trans- 
form is denned as p times the integral in (i.i), an additional 
term appears arising from this extra factor, p, and the simplicity 
of (2.8) is lost. 
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using (1.6), it appears that the essential content of Theorei"* 1 
II is simply the formulae (3.5) . . . (27) for putting? 
f(p}(g(p) into partial fractions in the case in which g(p) has 
no repeated factors. This is a well-known result,* fon~t 
for completeness we give a proof below. 

All of (2.5) to (2.10) are useful ; (2.6) is easily renaen*- 
bered by noticing that the coefficient of i/(p a r ) is tn^ 
result of putting p a r in f(p)lg(p), omitting the facto** 
(p a r ) of g(p] ; (2.7) is simply another way of writing? 
this statement. Similar remarks apply to (2.9) and (3.10)- 
The partial fraction formulae (2.5) and (2.6) are themselves 
very useful, cf. Example 4, below. With regard to (^.8) 
and (2.9), it is usually better to use (2.9) when g(p) has & 
few simple factors, as in the examples of this section - 
It sometimes happens that g(p) has complicated factors, 
which need not even be set out explicitly, as in the problems 
of 1 6, and then (2.8) is the more useful. 

To derive (2.6) we put y(p) into partial fractions by 
methods discussed earlier, by assuming 



Then 



Putting p = a r in this gives 

f(a r }-=^(a r -a^ . . . (^-^-O^-flH-i) (r n) 

r = i, 2, . . ., -, 

and substituting this value in (2.1 1) gives (2.6). To deduce 
(2.5), we notice that differentiating (2.4) with respect to 
p by the ordinary rule for differentiating a product gives 



Putting p = a r in this, gives 

g'(a r }=(a r -a^ . . . (a r -0 r -i)K-r+i) - (r-^)- (2. 

Using (2.12) in (2.6) gives (2.5). 

* Cf. Gibson, Treatise on the Calculus, edn. a (1906), 120. 
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Before proceeding further it is desirable to settle some 
questions of terminology. Firstly, suppose y(p) is given 
in the form (2.1). The simplest case is that of (24) in which 
g(p) is a product of n different factors, 

The numbers a : , 2 , . . ., a n are the roots of the equation 
g(P} - They may also be spoken of as the zeros of 
g(p), since they are the values of p for which g(p) is zero. 
Finally, they are the values of p for which the function 
y(p] is infinite, and in this sense are described as the poles 
of y(p). In the more general case in which g(p] has 
factors such as (p a^ n , a : is called a multiple (or more 
precisely an w-ple) root of g(p) = o, an n-ple zero of 
g(p\ or a pole of order n of 



Secondly, there is the ques- 
tion of the notation for a 
complex number. If 

z = x -f- iy 

is a complex number, we 
call x the real part of z, y 
the imaginary part of z, and 
x iy the conjugate of z. _ 
The number may be repre- < 
sen ted on the Argand diagram 
of Fig. i, and may also be 
specified by its modulus, \s\ \X(# 2 
(f> defined * by any two of 




FIG. i. 



the angle 



sn = 



V 



cos 



SC *V 

= -' tan ^ = 



We shall call this angle <j> the argument of x, written 
arg z ; other names for it are the amplitude of s, (amp #), 
or the phase of z. 

*'This definition still leaves <j) undetermined by 'an integral 
multiple of z-rr ; this does not matter in the applications made 
here. If desired, <f> can be restricted to the range it > < > IT, 
when it is called the Principal Value of arg x. 
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With the notation. (2. 1 3), any complex number x -\- 
may be written in the form 

ar+^V^+yy*, (2. 1 4-) 
where < = arg (x + iy). 

The importance of this in the present connection is 
when y(t) is found from y(p) by the use of (2.9), it 
frequently involve products or quotients of co 
numbers, and the most rapid way of reducing suclrx 
expressions .to a real form is by using (2.14) for each o 
the complex numbers (cf. Example 5, below). 

One pitfall of some importance must be mentioned 
the definition of </> in (2-13) is equivalent to * 

tft = arg( + iy) = tanr 1 ^/*), if x > o\ /. T ,- \ 
= TT + tan-^/a), if * < oj , ^ 5 ^ 

thus, if we had written <j> = tan~ 1 (y/w) in (2-14) (and this 
is very frequently done), there might be an error of 
in the argument of as.-}- iy, leading to an error in sign 
the final result. This point is particularly important 
Tables of transforms, and in problems such as Example 
below, where numerical values of some parameters 
not specified and thus may be taken with either sign. 
We now give some examples on the use of Theorem II- 

Ex AMPLE 3. V = ; ; TT ; TT -. . 

6 y Xp + iXP+ *)(P + 3) 
By (2.6), 



< 



and 3; follows from (1.6). Or, using (2.9), 
= - l ? -t:.. - 1 



* tan a is defined as the angle between Jir and TT whose 
tangent is x. 



II 

EXAMPLE 4. y = a , a . , *. 

' 2(^2 _|_ !)(2 + 4) 

Using (2.6) on y, regarded as a function ofp z , we find 

i . i 



Thus y = -t sin t -\ sin zt. 

4 3 24 

Problems on resistanceless circuits often lead to transforms 
which can be conveniently treated, as above, as functions 
ofp*. 

EXAMPLE 5. 

p -\- a 

y 1 

(P ~\~ fy[(P H~ j 1 *) 2 ~l~ ^ 2 ] 

^ (P + b)(p + n- in)(p +/* + in). (*' l6) 

This is a generalization of Example i, and could be 
dealt with either by resolving into partial fractions with 
real denominators, as was done in Example i, or by re- 
solving into partial fractions with linear complex denomin- 
ators, that is, by using one of the forms of Theorem II. 
Thus applying (2.9) to (2. 16) gives * 



y- 

y 



- by + 

-|- Conjugate) (2. 17) 
j & j \ /y 



Then writing 



* A pair of conjugate roots of g(p} = o give in (2.9) the sum of a 
pair of conjugate terms, of which only one need be written down. 
The method used below for reducing such a pair of terms to real 
form will be found very convenient. 
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and using (2.14) in the complex numbers in (2.17), we 
find 



wV re*(*+* l -#.) -k 

B r e i - = -- ^ Conjugate V 

2 J I 2z J J 



(a - b) 

r * _ 7 



.. 

/> 01 



Theorem II covers the case * in which g(p} has no repeated 
factors ; for example, it is not applicable to Examples 2 
or 4 above. It is easy to generalize the theorem to the 
case of repeated factors. Theorem II may be stated in 
the form that, to each linear factor (p <z r ) of the denomin" 
ator ofy(p) there corresponds a term 

fgf* . . . (MS) 

in the solution. 

The generalization is that, to each squared factor 
(p Z>) 2 of the denominator ofy(p) there corresponds a term 



( 

L 



i^e solution. 

* In (3. i ) it was stated that/() and g(p) were to have no common 
factors : this is desirable but not essential. With complicated 
expressions (such as those of 16) it would greatly increase the 
labour to separate out such a factor, and, in fact, if f(p) and g(p) 
both _have the same factor (p a), the term in e at corresponding 
to this in Theorem II has zero coefficient, and so the application 
of Theorem II gives the correct result. If g(p) contained a 
factor (p a) 2 , and f(p} a factor (p a), Theorem II could 
be applied after cancelling the factor (p a), whereas a cursory 
inspection would suggest that it is not applicable and that (2.18) 
would have to be used (this, also, would give the correct result). 
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And to each m-ple factor (p c) m of the denominator 
ofy(p) there corresponds a term 



in #ze solution. 

As an example, consider the problem of Example 2 
above, in which 

__ _ i _ _ _ i _ 
y ~ (P + *)*(P Z + 4)~(P + *Y(P - 2i)(p + aty 
Then by (2.19) and (2.18) 



+ 



leading to the same result as before. 

Finally, the question of " uniqueness " must be dis- 
cussed. To a given function y(t) there corresponds a 
single transform y(p) defined by (i.i). It is not clear, 
however, that to a given transform y(p) there corresponds 
only a single function y(t). Given y(p) we have sometimes 
been able to find a function y(t) which has y(p) for transform 
by looking up y(p) in the Table of transforms ; but there 
is so far no reason why the function of t found this way 
should be the only function which has y(p) for transform. 
The same question arises in the integral calculus when 
we seek a function which has a given function for differential 
coefficient, and there it is found that there are infinitely 
'many such functions (differing by additive constants). 

In the present case the matter is settled definitely by 
Lerch's Theorem* which states that, if two continuous 

* For a proof see C. and J., Appendix I. The actual theorem 
is a great deal wider than that stated above, and includes also 
the useful result that if two functions y^t") and y z (f) are continuous 
except at a number of ordinary discontinuities (i.e. jumps in value), 
and have the same Laplace transform, they must have the same 
points of discontinuity and can differ only at these points. 
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functions have the same Laplace transform they must be 
identical. Thus, if we find a continuous function y(t) 
from a given y(p), for example by the processes described 
above which use the Table of transforms, this is the only 
continuous function which has y(p) for transform. 

It follows also from Lerch's Theorem that, if k is a 
constant, and 

y(P) = > p > k, then y(t) = o, t> o . ( 2 . 2r ) 
provided y(t] is continuous ; and if y(t) is not continuous 
it can only differ from zero at isolated points. 

3. Some additional theorems 

THEOREM III. // y(p) is the Laplace transform * of 

y(f), and y(t] is continuous f and tends to y as t-> o, then 

\ the Laplace transform of dyjdt is py(p) y . (3.1) 

If, in addition, dyfdt is continuous and tends to y : as t-+ o> 
then 

the Laplace transform of -^ is p z y(p\ -py - ^. ( 3>2 ) 

Andifd s y/dt 3 . . . d n y/dt n are continuous and dy/dF -> y r> 
r ~ z >3> -, n ; i, as t-^- o, then 

the Laplace transform of -- 
J J dt n 



P n -y - ... -py n ^ - y n _,. (3.3) 

These results follow immediately by integration by 
parts. To prove (3.1) we have 



f .dy r -100 
o 6 dt dt = 



t ? the - functions concerned have Laplace 
f* 10 of , such as ', for which the 
esu t T c01 ? Verge h wever lar S e # is ch ^n, and 

not likS to , Ve - L t pla e transform s. Such functions are 
fs us 6 ^ pr blem to which the 

1 



ra ^' is> ]"* jUm P S in value > at 

extra terms arise corresponding to these points. 
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since ye~ Fi ->y as t-^-o, and ye~ vt ~^-o as ;- oo, 
provided p is sufficiently large. 
Similarly, to prove (3.2), we have 



Repeating the process w times gives (3.3). 

THEOREM IV. If y(p) is the Laplace transform of 
y(t), then 

f f i 

the transform of I y(t'}dt' is -y(p). (3-4) 

Jo P 

For 



roo ft r T ft -too T r 

e-^dt \y(f)dt'= \--e-v* \ y(t')dt' \ +- <r**y(t)dt 

Jo Jo" L p Jo Jo PJO 



since the integrated part vanishes at both limits. 

The results of Theorems III and IV can often be used 
to obtain new transforms : as a simple example it follows 
from (1.7) and Theorem IV that 

a f* . i 

13 - 57 is the transform of 1 sin at dt = ~(i cos at). 
XP +*) Jo ^ 

THEOREM V. // y(p) is the Laplace transform of y(t), 
and a) is a positive constant, then 

y ( ) is the transform of y(ajt). . (3.5) 

For f " e-**y(ti>t)dt = - 
Jo w 



This simple result is useful when it is desired to arrange 
the algebra of a problem in dimensionless form. For an 
example of its use see 14. 

4. The solution of ordinary linear differential equations 
with constant coefficients. 



l6 THE LAPLACE TRANSFORMATION 

Writing D n y for d n yjdt n , the problem is to solve the 
differential equation 



(4.1) 

where c lt . . ., c n are constants, /() is a given function 
of if, and j, Dy, . . ., D"" 1 ^ are to take the given values 
yw y - -, y n ~i when t = o. 

We multiply (4.1) by e~^ and integrate with respect to 
t from o to oo. This gives 

J V*{Dy + cj)-*y + . . . + c n y} dt = /(#), (4.2) 

where /(j>) is the Laplace transform of /(*), which can be 
written down* from Table I and Theorem I for the 
common types of function which usually occur. 

On the left-hand side of (4.2) we use Theorem III, 
which gives ' 



Or 



i}- (4-3) 

The equation (4.3) is called the subsidiary equation 
corresponding to the differential equation (4.1) and the 
given initial conditions. It is formed from the differential 
equation (4.1) by replacing D r y, r = o, ...,, on the 
left-hand side, by fry- on the right-hand side, f(t) is 
replaced by its transform f(p), and terms containing the 
initial conditions are added according to the rule : 

for Dy we add y 

for D*y py + yi 

foTWy ptyo+pfr+yt 

for D^ . p'-iy, + y-^ + ....+ y^. 

*The solution for a function whose Laplace transform is 
not known, or even which has no Laplace transform, can be 
obtained by the use of Theorem IX, 24. 
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The subsidiary equation (4.3) gives y(p), the transform 
of the solution y(t) of the problem. To find y(t) from 
y(p), the methods of 2 are used. Since, by the unique- 
ness theorem of 2, only one continuous function y(t) 
corresponds to a given y(p), it follows that this will be the 
unique solution of the problem. 

EXAMPLE i. (D 2 + 1)3; = o, > o 
with y = y Q , Dy = y ly when t = o. 

The subsidiary equation is 



Thus 



And y = y Q cos t -f y\ sin f. 

EXAMPLE 2. (D 2 + 4D -f B)y = i, * > o 
with 3> = o, Dy = i, when t = o. 

Since the transform of i is ifp, the subsidiary equation is 



The function y corresponding to this has been found in 
2, Ex. i, and for a similar but more general case in 
2, Ex. 5. 

EXAMPLE 3. (D + i) 2 y sin zt, t > o, 
with y and Dy zero when t = o. 

The subsidiary equation is 



The function 3; corresponding to this y has been found in 
2, Ex. 2. 

EXAMPLE 4. (D 2 + n z )y = cos /, i > o, 
with .y == 3/0, D^v == ^ 1; when t o. 

2 
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The subsidiary equation is 



Thus y = _ t _ + 

y 2 +7Z 2 ) 2 ^ + W . 

Using (i.n), (1.7), (1.8), it follows that 

t . v, 

y = sm ;rf + 3> cos nf + sin 



EXAMPLE 5. (D + i)y = t*e-*, t > o, 
with y = y 0> w hen * = o. 

Using (1,5) and Theorem I, the transform of t z e- t 
is found to be 2/(p + i) 3 . Thus the subsidiary equation is 



2 

(7TT? 






2 , y* 



. 
Then, using (1.5) and Theorem I, we find 

y = $t*e-*+ yo p-4. 

5. Simultaneous ordinary linear differential equations with 
constant coefficients. 

Here, proceeding as in 4, using the transformation 
procedure on each of the equations in turn, gives a 
system of algebraic equations, the subsidiary equations, 
for the transforms of the solutions. Suppose, for example, 
it is required to solve a system of n second order equations 



r= i, 2, . . ., n, t> o, (5.1) 
where the /,(*) are given functions of t, and y r and Dy r , 

T 7/1 ' ' '' W> are to take the values u r and w r when = o. 
Multiplying the equations (5.1) by -*, integrating 
with respect to t from o to oo, and using (3.1) and (7.2) 
gives the system of n subsidiary equations 
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in 



+ K(M + w ) + M S ], r = i, ...,. (5.2) 
j-i 

These have to be solved for y lt . . ., <y n) an d then 
y lt . . ., y n are determined as in 2. 

EXAMPLE. (D 2 + z)x Dj/ = 
D^-h(D 2 + 2)3; = 

to be solved with * = tf , D = y = Dy = o , when i = o. 
The subsidiary equations are 



px + (p 2 + z)y = ! . 
Solving we find 



I 2X I 



4) 3 V + i 2 + 4/ 

Thus x = 4- &(2tf i) cos ;? + ^(2 i) cos zt 

6. Comparison with the ordinary method of solving differential 
equations, and with Heaviside's operational calculus. 
The ordinary method of solving the linear differential 

equation of order n, 

<j>(D}y = (D n + cjD"" 1 + . . . + c n )y =/(*)> (^ I ) 
consists of ^finding the complementary function 

4**', . . (6.2) 

I r-l 

where Aj, . . ., A n are arbitrary constants, and the <%, 
. . ., a n are the roots of <j>(p) = o (provided these are all 
different, if they are not the complementary function 
takes a slightly different form provided for by the theory) 
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and adding to this a particular integral determined from 
</>(D) and /() by certain rules. This gives the general 
solution of (6.1). To find the solution satisfying given 
initial conditions, it is necessary to solve n algebraic 
equations for the constants A 1? . . ., A n . 

In the Laplace transformation method, the subsidiary 
equation (4.3) is 

y = T7J(f(p} + terms involving the initial conditions}. (6.3) 

To determine^ from (6.3) we have to put the right-band 
side into partial fractions,* and for this the roots a^ . . ., 
a n of <j)(p) = o are needed ; thus these must be calculated 
in either method, but when they have been found, the 
work by the Laplace transformation method is much 
shorter ; for example, if a^ . . ., a n are all different, 
y can be -written down by (2.8), while with the older method 
of solution the n algebraic equations for the constants 
Aj, . . ., A n of (6.2) still have to be solved. Thus the 
advantage of the Laplace transform method over the 
ordinary one increases enormously as the order of the 
equation to be solved increases. 

The same remark applies to simultaneous linear differ- 
ential equations the greater the number and order of 
these, the greater is the advantage of the Laplace trans- 
formation method. _ 

The part of y(t] derived from the term /(p) /</>(/>) of 
(6.3) is a particular integral of (6.1), but it may differ from 
the particular integral found by the rules of the older 
theory by terms of the complementary function. 

To see the connection between the present method 
and Heaviside's operational method, we again discuss 
(6.1). Heaviside usually considered only the case in which 
the initial conditions were zero, and f(t) was his " unit 
function ", which is zero for t < o and unity for t > o, 

* If /CO is not one of the simple functions dealt with hitherto, 
it may not be possible to put /()/<() into partial fractions, but 
the solution can always be found as in 24. 
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and is represented by the symbol 1. On the left-hand side 
of (6.1), D is replaced by p ; this is regarded as an operator 
(not as a number as in the Laplace transformation method). 
Thus (6.1) leads to 



t ' ' ' (6<4) 

which is the operational expression for the solution y. 
The solution itself is then derived from (6.4) by Heaviside's 
expansion theorem which is similar in type to (2.8). 
For this problem, the Laplace transform of the solution is 

(6 ' s) 



and this differs from (6.4) only by the extra factor p in 
the denominator. In all cases the Heaviside operational 
expression is just p times the corresponding Laplace 
transform. 

As remarked in i, the Laplace transformation pro- 
cedure is designed for, and applies immediately to, the 
" initial value problems " of dynamics and electric circuit 
theory, in which differential equations have to be solved 
for times t > o with given initial conditions at t = o. 
In other fields, such as the theory of elasticity, " boundary 
value problems " arise, in which a differential equation 
has to be solved in a region, say o < x < /, and the 
solution has to satisfy some conditions at the boundary 
x = o, and some at the boundary x L The older 
method of solution, as the sum of a complementary 
function (6.2) and a particular integral, applies equally 
well to both types, since, in a boundary value problem, 
there will still be n conditions, some at x = o, and some 
at x I, to determine the A 15 . . ., A n . The Laplace 
transformation can still be used, since the initial values 
y lt . . ., y n may be regarded as arbitrary constants ; some 
of these may be fixed by the conditions at x = o, and the 
remainder can be determined from the conditions at x = / 
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As an example, we consider the deflection * of a uniform, 

uniformly loaded, beam of length Z, freely hinged at the 

origin, and clamped horizontally at the other end, x== L 

, Writing D for differentiation with respect to x, the differ- 

ential equation for the deflection is 

EID 4 jy = to, o < x < I, (6.6) 

where E, I, and w are constants. This has to be solved 
with x = DV= o at x = o, and x = D# = o at x I. 
We solve (6.6) with values o, y^, o, y 3 of y, Dy, D z y, D 3 ;y, 
respectively, at x=.o, regarding y^ and y 3 as unknown 
constants to be determined from the conditions at x I- 
Writing y for the Laplace transform of y with respect 
to the present independent variable x, the subsidiary 
equation is 



mi 

Thus 



The conditions y = Dy = o at x = I, give 



24EI ' 6 



mr 

Thus 



and y = 

7. TAe justification of the solutions of 4, 5 /row 

pure-mathematical point of view. 

In deriving the subsidiary equations of 4 and 5 some 
assumptions were made, which, from the point of view 

* This method is well suited to problems of this type, par- 
ticularly with concentrated or variable loads. For other examples 
see 21 ; Jaeger, Math. Gazette, 23 (1939), 62 ; Pipes, Journ 
App. Phys., 14 (1943), 486. 
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of strict pure mathematics, would have to be defended. 
For example, it was assumed that y and its derivatives 
had Laplace transforms. To make the solution quite 
rigorous further discussion would be needed, but this 
can be avoided by checking the results obtained by the 
methods above, that is to say verifying that they do satisfy 
the differential equations and initial conditions of their 
problems. This can be done quite generally, and since 
the process itself is purely algebraical and adds nothing 
to the theory of the subject it is sufficient to state the 
results here. 

For the problem of 4, the ordinary linear differential 
equation of order n with constant coefficients, it can be 
verified * that the solution obtained by the Laplace 
transformation method does satisfy the differential equation 
and initial conditions. 

For the problem of 5, simultaneous ordinary linear 
differential equations, it is possible to give a similar general 
verification that the solutions obtained by the Laplace 
transformation method satisfy the differential equations, 
and that they satisfy the initial conditions, provided that 
the determinant of the coefficients of the highest powers of 
D in the terms does not vanish. This restriction is re- 
quired, not because of a flaw in the method of solution, 
but because of a fault in the specification of the problem. 
To see what it implies, consider the system of equations 
(5.1) : in this case the determinant in question is 



If this vanishes it is possible, by adding together appropriate 
constant multiples of, the equations (5.1), to form at least 
one equation containing no second derivatives. Thus, 
instead of having a system of n second order equations, 
which could be solved with 2n initial conditions, we have 
in fact a system of n i (or less) second order equations, 
* C. and J., 34, 35 ; Doetsch, loc. cit., Kap. 18. 
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and some first order or algebraic equations, and for such. 
a system only zn i, or less, arbitrary initial conditions 
can be prescribed. If the problem is stated with 
initial conditions, there must be relations between these 
or there will not be a solution. 

Problems of this type occasionally arise in practical 
applications, also similar difficulties occur in switching 
problems in circuit theory. These questions are discussed 
in 20. 

EXAMPLES ON CHAPTER I 
i. Verify the following Laplace transforms : 



m 


,w 


i 


|V<T* . /i*^ 


p 


(a - *>r ; 


(p- a)(p- Z>) 
i 


(a ft) 
. (i cos at) 


p(p Z + <2 Z ) 
I 


I 




P 


ab(b z a 2 ) 


(p z + a i ~){p' 2 ' ~\- b z } 
i 




Pl(P + Y + W 


" + ' s ^;:^ + , ra 



z. Solve the following differential equations with initial values 
3>o> yi, - , of y, Dy, etc. 

(i) (D + M )y = i. fy = -+ fyo- 

L, f* \ A 1 

(ii) (D 2 + 5 D + 6)3; = i. 
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(iii) (D 2 + n*)y - sin o>f. 

-M 2 ) . . i . ."1 

Sm nt + y C S " - ttn ' ' 



(iv) (D + i) 2 y = sin f. 



3- If d*yldx* = x, o < x < /, 

with 3; = dy/dw = o at = o and x I, 

show that y = 2 ( 3 + zl a ~ sZ 2 ). 



28 THE LAPLACE TRANSFORMATION 

\ 

series with a leaky condenser of capacity C and leakage 
conductance i/G, the subsidiary equations are 



Returning to (8.6), we remark that it involves the 
initial current in the inductance and the initial charge 
on the condenser (these are related to the initial kinetic 
and potential energies, respectively, of the current). 
Further, since the terms involving them appear in the 
right-hand side added to the transform V of the applied 
voltage, it follows that the effect of initial charge Q on 
the condenser is the same as that of an applied voltage 
whose transform is ( Q/Cp), that is of a constant voltage 
( Q/C) ; similarly, the effect of initial current 1 in the 
inductance is the same as that of an applied voltage whose 
transform is Ll, this will be found in 19 to be an im- 
pulsive voltage LlS(f). The same remark is true in. 
general that initial conditions may be treated as voltage 
sources of the appropriate types. 

9. Examples on L, R, C circuits. 

EXAMPLE i. Constant voltage, E, applied at t o to 
L t R, C in series,* with zero initial current and charge. 

Here the subsidiary equation (8.6) is 



Thus, 

T f \ 

= * + (R/L)p + (i/LC)) = ' ' ^ 



R 2 _ * 2 _ J 

~LC ^~LC 
The notation (9.3) will frequently be used later. 

* For the exceptional case L = R = o and related problems, 
see 20. 
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From (9.2), using Theorem I and (1.7), (1.5), and 
(1.9), respectively, we find 
p 

I = -^-e-P* sin nt, if n 2 > o 
nLi 

TT 

= r flr' 1 *, if n 2 = o 

JL 

E 
= Tf-e-M sinh fe, if n 2 = k z < o. 

RL, 

EXAMPLE 2. .4 condenser of capacity C charged to 
voltage E is discharged at t = o through an inductive re- 
sistance L, R. 

Here Q = CE, I = o, V = o, and (8.6) becomes 



which is the same as (9.1) except for the negative sign. 
This implies that the current starts flowing away from 
the positive side of the condenser. 

Using (8.5) and (9.4), the transform of the charge on 
the condenser is found to be 

CE CE ' 



where the notation (9.3) has been used. Then, in the 
case n z > o, we have 

Q.= CEfi-Wcos nt 4- ^ sin nt\. 

EXAMPLE 3. Steady current E/R is flowing from a battery 
of voltage, Ej through a parallel combination of condenser 
C and inductive resistance L, R (Fig. 3), when at t o the 
switch S is opened. It is required to find the potential across 
the condenser terminals. 

After the switch S has been opened, the circuit consists 
of L, R, C in series. The initial charge Q on the condenser 
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is CE, and the initial current 1 in the inductance is 

the minus sign being chosen, since, from the conditions 

of the problem, this is flowing away from the P os1 ^ 

side of the condenser. With these initial values (8.0) 

becomes 



1 = 



where the notation (9.3) has been used. Hence, con- 
sidering only the case n 2 > o, 



nt + - sin nt\ . 




FIG. 3. 

If v = Q/C is the voltage across the condenser ter 
minals, it follows from (8.5) that 



__Q_E_ 

v c~p 



R) 



+ 



+ i) 



cos nt 



-( -= -rr^ 
n\zL RC 



(p + ju) 2 + a 



sn 
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This is the simplest type of circuit breaker restriking 
problem in which it is required to find the voltage across 
the contacts of a switch after they have been opened. 
In this case the voltage across the switch S is E v. 
Anothe-r method of solving this problem is given in Ex. 2 
at the end of this Chapter. 

EXAMPLE 4. Alternating voltage E sin (o>t -f- <) applied 
at t = o to an L, R, C circuit with %ero initial current 
and charge. 

Since sin (wt + 0) = sin tot cos -f cos wt sin <^, its 
transform is 

p sin <^> + w cos <j> , * 

2 + o> 2 ' ' ' (9 ' 5) 

Thus, using the notation (9.3), the subsidiary equation 
(8.6) gives in this case 

. _ E j> 2 sin $ -{- tap cos <ft , ,. 

L ' ' * (9 ' 



We consider only the most interesting case w 2 > o. 
Here to find I we may either put (9.6) into partial fractions 
with real quadratic denominators, or notice that its 
denominator is 



(p i(jo)(p + iot)(p + M- in )(P 
and use (2.9). The latter method gives * 

, EC o) a sin<^ + ia> 2 cos<j> ico . . n . ^} 

I = H r ~fv i \8 t sr 1 * 2 + Conjugate V 
LI 3zo)[(ju, + W) 2 + w 2 ] j & j 

E - + m 2 sin <> + a>- + cos 



f 

""LI 



ut+tni 



Conjugatej (9.7) 



* Notice that if only the steady state part of I is required, only 
the terms of (2.9) arising from the zeros i<a of the denominator 
of (9.6) need be written down. The other zeros /u, in give 
transient terms. See also 13. 
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To reduce this to real form, -we use (2.14) as iti 2, 
Ex. 5. Thus, writing 

-[a/to) + i(u>* -V?- n z )] = Z", . (9- 8 ) 

O) 

(/i 2 + w 2 2 ) + atjun = Z^, 
[(/x 2 2 ) sin ^ |u,aj cos <j -+- i[no) cos < 2^n sin ^>] 

= 2>2^ *| 

(9.7) becomes 

S 2 ), (9.9) 



where Z and 8 may be written out explicitly as the mo d.xilus 
and argument of the complex number on the left of C9*^)i 
and similarly Z lt B lt Z 2 , 8 2 . 

10. Electrical networks 

It is convenient to regard a complicated electrical net- 
work as built up of a number of circuit elements of th< 

L 5 U 




M rs 



*^OT^ W 1 




FIG. 4. 

type studied in 8, with their terminals connected in som< 
specified manner. Let L r , R r , C r be the inductance 
resistance, and capacity in the r-th of these elements 
(Fig. 4), let V r be the voltage drop over it, o let I r be thf 
current in the element, with initial value I r , and let Q 
be the charge on the condenser, with initial value <!$,. 
To he quite general mutual inductance must be included 
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so let M Ta be the mutual inductance * between the r-th 
and s-th elements. 
Then the equations for the r-ih circuit element are 

JT rk _ j-r 

L. _> M ? = V do i.} 

r ^ iVJ. rs . v r , . ^IU.A; 



dt =I r . (10.2) 

Proceeding as in 8, the subsidiary equation for the 
y-th circuit element is found to be 

#1 4-2-s 1 1 ="V 4- Li 4- 5! M 1- 



where z r = L r + R r + --, . . (10.4) 

^"rP 

*rt = M ra p (10.5) 

From the form of the network, further equations can 
be written down connecting the 1 T and V r in different 
circuit elements. Firstly we have Kirchhoffs First 
Law f 

]^] I = o, at a junction. . . (10.6) 
This leads to 

^ I = o, at a junction. . . (10.7) 

Kirchhoff's Second Law states that, for any closed circuit 
in the network, the sum of the externally applied voltages, 
V, in the closed circuit, equals the sum of the voltage 
drops, V r , over the circuit elements included in the closed 
circuit. That is, writing ^ to denote summation round 

o 

* M r , will have a sign depending on the way in which the 
coils are wound, and on the convention of sign we set up for 
the Currents in the different circuit elements. 

t Usually labour can be saved by specifying the currents, 
either: directly or by the 'use of mesh currents, so that they satisfy 
(10.6) automatically. 

3 
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a closed circuit V y r _ V y, 

CO 

and thus ZV r = ]T?. ' ' 



Using (10.8) in conjunction with (10.3) we have 



By choosing suitable closed circuits, sufficient equations 
of types (10.9) and (10.7) can be found to give all the I r . 

When voltage V, applied' to the terminals of any net- 
work with zero initial conditions, produces at the terminals 
current I given by 

V, . . . (10.10) 



we call s the " generalized impedance " of the network 
viewed from these terminals ; * for example, % r defined 
in (10.4) is the generalized impedance of an ii r , R r> C r 
circuit. Now for a single circuit element of a network 
with no mutual inductance and with zero initial conditions, 
(10.3) becomes 

* r I r = V f . . . . (io.ii) 

That is, I r , V r , and z r are connected in the same way as 
current, voltage, and resistance in direct current theory. 
Thus in this case the generalized impedance of the network 
can be written down by combining the z r of the circuit 
elements which comprise it according to the ordinary laws 
for combining resistances in direct current theory. For 
example, the generalized impedance of a combination of 
n generalized impedances z 1} z 2 , . . ., z n in parallel is 
given by 

1 = 1 + 1 + ...+JL. . (10.13) 

^ 



* 



* 2(10)) is the ordinary complex impedance of alternating 
current theory. 
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Again, the generalized impedance of the ladder network, 
Fig. 5, can, by repeated application of the result for two 
resistances in parallel, be written down as the continued 
fraction 

*==*!+- - 1 - . . (10.13) 



, 









FIG. 5. 

It follows in the same way from the remarks at the 
end of 8, that the current in any branch of a network, 
due to initial current 1 in an inductance and initial charge Q 
on a condenser in that branch, is given by 

= Ll - (QfCp) . . (10.14) 



where z(p) is the generalized impedance of the network 
viewed from a pair of terminals in the branch. 

n. Examples on electrical networks 

EXAMPLE i. The two circuits of Fig. 6 are coupled by 
mutual inductance M. Voltage V(t) is applied to the 
primary at t = o with zero initial conditions. It is re- 
quired to find the secondary current (the primary resistance 
is neglected). 

Here the subsidiary equations (10.3) give 



Mpl, 4- L^ + R 2 + JLl a = o. 
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FIG. 6. 
If V = E, constant, so that V = E/p, this becomes 

ME 



ME 



and 

where a = 



sin 



a". 



The case L A L 2 = M 2 is discussed in 20. 

EXAMPLE 2. Ira Z/ze circuit of Fig. 7 Z/ze condensers are both 
charged to voltage E, when, at t = o, the switch S is closed. 

A I L B L 
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We choose the currents as shown in Fig. 7 to satisfy 
Kirchhoff's first law. Then the subsidiary equation 
(10.9) for the closed circuits ABD and BGHD are, 
respectively, 



CE CE 
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Solving for I, we find 
E ' Zj r 2n * 



T= 

~ 



where 2 = i/LC. 
Therefore 



E f i + VS 
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EXAMPLE 3. In the circuit of Fig. 8 it is required to find 
the current I in the inductance L caused by prescribed 
current I supplied at the terminals AB, starting at t = o 
with zero initial conditions. 

1 
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>R = 


= c | 
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A A I 


1 



FIG. 8. 



From (10.12) it follows that the generalized impedance 
of the combination of L, R, C in parallel is 



RLp 



* ^ R ^ Lp RLC 2 + Lp + R. 

Also the generalized impedance of the inductance L i? 
Lp. Since the same potential is applied to both we hav 
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For example if I = i, for i> o, so that I = i/p, 



and, if ^ > o, 



where a - 



p(RLCp z + Lp + R)' 
er** sin (/ft + ^ 



- JL _ a 2 , = tan" 1 (]8/a). 



EXAMPLE 4. Bridge networks. 

We consider a general bridge network represented by 
the scheme of Fig. 9, and suppose known voltage V to be 



-Is 




H 



FIG. 9. 



applied at i = o with zero initial conditions. It is required 
to find I 5 , which is the galvanometer current in most 
applications. Choosing the currents as shown in Fig. 9 
to satisfy KirchhofFs first law, the subsidiary equations 
for the closed circuits ABDHA, ABFA, BDFB, re- 
spectively, are 



= o 
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Solving for I 5 we find * 

1 +#2 O 
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where A = z^ 



# 



# # ^ 



The way in which the bridge can be used for measure- 
ments is determined by the form of z z s 3 si-jjs^ in (11.4) 
as a function of p. If, by suitable choice of the circuit 
constants, this can be made to vanish for all p, we have 
1 5 = o, and thus, by (2.21), 



whatever the applied voltage V. As an example of this 
case consider the bridge in which x consists of inductance 
L and resistance R x in series, z z of resistance R 2 , # 3 of 
resistance R 3> and # 4 of capacity C and resistance R 4 in 

parallel. Here 



*! = !# + R 15 # 2 = R 2 , % = R 3 , # 4 = 



R 4 Q>" 



TTniQ y y y y 

JL J..LU.O ***2 3 14 

Then if 
and 

we have I 5 = o for t > o, whatever the form of the applied 
voltage may be. 

* It is hardly necessary to point out the advantage of using 
determinants, particularly when, as here, only one of the currents 
is needed. 
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,As an example of a bridge whiqh cannot be balanced in 
this way, consider the case in which 

z 1 = ltp -\- RI, #2 K-2) #3 = RS> 

and # 4 consists of resistance ^ in series with a parallel 
combination of resistance r z and capacity C, so that 



and 2 2 #3 



_ 

r " Q>+1 (i i. 6) 

Clearly (11.6) cannot be made to vanish for all values of p 
by choice of the circuit constants. However, we rnay 
use the circuit as an alternating current bridge by requiring 
that the steady state part of I 5 is to vanish when the applied 
voltage V is sinusoidal and of frequency o>/27r. Antic- 
ipating the results of 13, this requires that the real and 
imaginary parts of the numerator of (n.6), with j& replaced 
by ioi, shall both vanish, that is 

Lfa + r a ) = 2 (R 2 R 3 - Vi) \ (T x~\ 

and LCr^co 2 = Rj^ + r a ) - R 2 R 3 J ' v . " x/ 

Alternatively the bridge may be used as a ballistic 
bridge ; in this case a battery of voltage E is applied at 
t o with zero initial conditions, so that V = E/p in 
(11.4), and the requirement is that the total charge passed 
through the galvanometer is to be zero, that is 

TOO 
Ldt = o. 



i 5' 
'0 



By Theorem VIII, 22, this requires 

lim 1 5 = o. . . . (n.8) 

p->0 ^ 

It is easy to verify from (11.5) that A does not vanish 
when p = o, so from (11.4) with V = E/p, and (u.6), 
it follows that for (n.8) to be satisfied we must have 
R 2 R 3 = R^ + r z Y 
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In practice, instead of applying a battery to an initially 
dead circuit, it is usual to open the battery circuit when 
steady current is flowing from the battery. In that case 
the impedances %, . . ., # 4 of Fig. 9 must be precisely 
specified, and the circuit of Fig. 9 (with the arm AHD 
removed so that I 2 = Ij) has to be studied with known 
initial values of the currents in the inductances and the 
charges on the condensers. 

12. The natural frequencies and damping factors of a circuit 
The transforms of the currents, say I l5 . . ., I n , in the 
various branches of an electrical network are found by the 
solution of equations of types (10.9) and (10.7), and 
always take the form 

t _ fr(P) _ T (,~ ,\ 

r ~~' ' ' ' " ' ' ' ( } 



where <3>(p) is the same for all l r , it is in fact the determinant 
of the coefficients of the l r in the equations and depends 
on the circuit only. The numerators f r (J>] depend on r, 
and involve the initial conditions and the applied voltages. 
Before we can find the I r we need the roots of the " period 
equation " 

o, . . . (12.2) 



and we know that to a root ft *" of (12.2) there cor- 
responds a term e~^ sin (vt -j- 8) in the I r , that is, a term 
of period 2w/v and damping factor er 1 **. 

The essential point is that the period equation (12.2), 
whose roots give the natural frequencies and damping 
factors of the circuit, is always found in the course of the 
calculation -of any of the l r , but the I r contain in addition 
the complete solution of the transient problem. Thus in 
ir, Ex. 2, it follows from the first expression for 1 that 
the natural frequencies of the circuit are n[(3 Vs)]*/ 377 - 

In demonstration problems, such as most of the examples 
solved here, 3>(p) is usually arranged to have obvious linear 
or quadratic factors. In practice $(p) may be a poly- 
nomial in p of degree n with coefficients given numerically, 
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an example of such a problem is given in 14. Here 
we give some references on the methods available for the 
solution of algebraic equations. 

For cubic and quartic equations algebraic methods are 
available, described in text books on algebra.* For the 
cubic, some tables and charts are given in Jahnke-Emde.'j 1 
Numerical methods can be used for cubic and quartic 
equations, and must be used for equations of higher degree. 
Real roots can be found from first principles, using a graph 
and improving the rough values obtained in this way, 
either by interpolation (an example of this is given in 14) 
or by Newton's method (W. & R., 44). A cubic has 
always one real root which can be found in this way ; 
dividing by the corresponding factor leaves a quadratic. 
The more sophisticated methods for finding real roots are 
Homer's -method (C. S., 475, W..and R., 53) and the 
root squaring method (W. & R., 54), the latter is prefer- 
able. Problems on resistanceless circuits lead to equations 
in p z with real (negative) roots which are easily solved by 
the above methods. For equations having several pairs 
of complex roots, the modified root squaring method seems 
the most satisfactory ; J another method, and also a method 
of successive approximation to a complex root, is given, by 
Frazer and Duncan. 

If any of the roots of the period equation have positive 
real parts, the solution will contain terms which increase 
exponentially with the tune ; in this case the system is 
said to be unstable. If all the roots have negative real 
parts the system is stable. Clearly it will be useful to 
have a result which will indicate whether a system is stable 

* Charles Smith, Treatise on Algebra, Edn. 5 (1921), referred 
to as C. S. ; Turnbull, Theory of .Equations (Oliver and Boyd, 
1939)- Whittaker and Robinson, Calculus of Observations 
(Blackie, 1924), referred to as W. & R. 

t Tables of Functions (Edn. 2, Teubner, 1933). 

I Brodetskyand Smeal.Proc. Cambridge Phil. Soc., 22 (1934), 83. 
Note that when there are several pairs of complex roots the method 
given in W. & R. is unsatisfactory. 

Proc. Roy. Soc., A 125 (1929), 68. 
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or unstable, without the labour of solving the period equa- 
tion. One such criterion is Hurwitz's * which states that 
the roots of the equation 

all have negative real parts if all the determinants 



a* a* 



a, a* a* 



, . . . 


% o 
3 a 2 


... 

flj O.Q O 






iZo__1 dn 








an JL 


(12 


3) 



are positive, where, in writing down the determinants, 
a r is put equal to zero if r > n. 

13. Steady state alternating currents 

It was remarked in 12 that if voltage V is applied to 
any circuit at t o, the transform of the current I r in any 
branch of the circuit will always take the* form 



(13-1) 



where f r (p) and $(_) depend on the circuit only, and 
g r (p) involves the initial conditions. Now suppose 

V = E sin (cat + </>,) 

cos d> -\- p sin d> , x 



so that V = 

and (13.1) becomes 



cos ^ + p sin <f>)f r ( p) _j_ gr(p) 



, 

' JJ 



If the system is a dissipative one, the roots of <(p) = o 
will all have negative real parts, and thus will give rise to 
transient teams. If only the steady state altern ating current 

* Frank- v. Mises, Differentialgleichungen der Physik (Edn. 2, 
1930), Vol. i, p. i6a. This is equivalent to Routh's Rule, 
Advanced Rigid Dynamics (Macmillan, Edn. 4, 1884), Chap. vi. 
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is required, we need only evaluate the two terms of (^ - 9) 
corresponding to the zeros i *<*> f P* + Ctj2 - These 



2icu$(i(a) 

' E 



Conjugate! 

J 



where 



^Thus, if the transform of the current in any part o a 
circuit due to applied voltage V and any initial conditions 
has been found, the steady state alternating current can 
be deduced immediately. An example has been given in 
(11.7). As another we find the steady state alternating 
current in the secondary of Fig. 6 due to the primary 
voltage E sin (cot + 0). Here (n.i) and (13.3) give 
Mj>E (a) cos <f> + p sin 



_ 

(p* + o^KI^L, - My + R 2 L^ + 

The steady state alternating part of I 2 is , 

' Mi'a>E(o) cos ft + *'a> sin </>)g icof , ^ . \ 

J 



/ . , . SN 

COS (ait + S), 

where Z = {R a V a +[(L a /C 2 ) - (L a L 2 - M 2 )^ 2 ] 2 }* 



and 8 = arg {[(L^C,) - w^L, - M 2 )] + z 

^ If 5f(p) is the generalized impedance of a dissipative 
circuit looking in at a pair of terminals, and voltage Ee* 60 * 
is applied at these terminals with zero initial conditions, 
the transform of the input current will be 

E 

(p-io>}z(p) y ' 
and the steady state current will be 
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14. Numerical treatment of more complicated circuits 

In the problems previously considered, the denominators 
of the transforms of the solutions have had obvious linear 
or quadratic factors, so that explicit expressions for the 
solutions could be written down immediately. In more 
complicated circuits, the denominators may contain ex- 
pressions of the third or higher degree in p which have 
no simple factors. If the circuit constants are given. 
numerically, the zeros of the denominator can be found 
by one of the methods mentioned in ra, and the solution 
completed ; if the circuit constants are unspecified, families 
of curves can be obtained which will show the behaviour 
of the circuit under any circumstances. 



+L 



FIG. 10. 



In this section we study, from both these points of 
view, the current I in the circuit of Fig 10, due to unit 
volTage applied at t = o with zero mitial conditions. In 
the usual way we find 

(14-1) 




As an example in which th ^ c = 2 , x 10 - 

;*Xi *='.o- ohms - Then 

becomes - 4 O 2 + i 10 ) 

To evalulte l! tCoots ^t^^^. 



.1) 
(14.2) 
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are needed. Putting p io 5 # in this, gives the more 
manageable equation 

/(a) = x 3 -f i-4* 2 + x + 0-4 = o. . (14-4) 

This cubic can be solved by any of the methods of jj 12, 
but it is very simple to work entirely from first principles. 
From a graph it is found that there is a real root at 
approximately x 0-77. This root can be found 
more precisely as follows : we find /( 0-77) = 0-00353, 
/( 0-78) = 0-00279, so tnat tne root ^ es between 
0-77 and 0-78 ; further, if /(#) ran linearly between 
its values at these two points it would vanish at 0-776- 
By repeating this process, as many places of decimals as 
desired can be found. Working with three places of 
decimals only, we divide f(x) by (x -f 0-776), and obtain 

f(x) (x -f o'776)(;x; 2 + 0-624*: + 0-516) 

= (* + o- 77 6)[(* + o- 3 i2) 3 + (o-6 47 ) a ]. (i4'5) 

Returning to (14-2) we find, u'sing Theorem II as in 
2, Ex. 5, that the function whose Laplace transform is 



s 



where 8= a r g (P+a-7j 2 +2^)+arg(a-^+^). (14.7) 



Using the values a = 0-776 X io 5 , ^ = 0-312 X io 5 , 
V) 0-647 X io 5 , a = io 10 , found above, and writing 
t io- 6 T, so that T is the time in micro-seconds, (14.2) 
and (14.6) give finally 



sin ( . 6 47 T- 1-485) 

(14.8) 

To plot I from (14.8) we notice that o-o647T = 77/8 
when T = 6-07, so that steps of 6-07 in T correspond to 
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increases of 22^ in the argument of the sine, and only 
eight values of the sine need to be looked up. The cal- 
culation may be set out as follows, A and B being the first 
and second terms on the right hand side of (14.8), 



T 


o-o776T 


*- 


A 


0-03I2T 


f&' 


r)t-d' 


sin 


B 


10*1 





o 


I 


2-527 





I 


-85-1 


0-996 


-1-527 


i-ooo 


6-o7 


o-47i 


0-624 


1-577 


0-189 


0-828 


62-6 


-0-888 


1-127 


0-450 


12-14 


0-942 


0-390 


0-986 


0-378 


0-685 


-40-1 


-0-644 


0-676 


0-310 


18-21 


1-413 


0-243 


0-614 


0-567 


0-567 


-17-6 


0-302 


0-263 


o-35i 


24-28 


1-884 


0-152 


0-384 


0-756 


0-470 


4'9 


0-085 


0-061 


0-445 


30-35 


2-355 


0-095 


0-240 


0-945 


0-389 


27-4 


0-460 


0-274 


0-5^4 


36-42 


2-836 


0-059 


0-149 


1-134 


0-322 


49-8 


0-765 


0-378 


0-5*7 



It very frequently happens in practice that it is not the 
transient behaviour of a circuit for given numerical values 
of the circuit constants which is needed, but a general 
description of the transient behaviour for any values of 
the circuit constants. The circuit of Fig. 10 will now 
be studied from this point of view. The first thing to 
do is to express (14.1) in a form involving only simple 
and fundamental quantities, such as dimensionless ratios, 
time constants of parts of the circuit, etc. Thus, writing 



(IA. 1} becomes 
V 



(i4-io) 



^-pt+dp*.. . _ 
A further simplification comes from expressing (14.1) u 
terms of A = p/i, so that the solution is obtained in terms 
of nit by Theorem V. Thus writing 

k = nln i , . . (14-") 



+ A 



i?T 
we have jRI = ^ 

The right hand side of (14.12) contains the two dimen- 
sionless, parameters b and k ; to study the general behaviour 
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of the circuit we can regard one of these as fixed at some 
arbitrary value, and see how the solution changes 'as the 
other takes all possible values. We shall give b fixed values, 
and for each of these allow k to run from o to oo, that is, 
say, R to decrease from oo to zero. 

First we have to study the position of the roots of the 
period equation 

x + k = o . . (14-13) 



for constant b and k >o. We could give k a number of 
values and solve (14.13) for each of these, but this would 
be very laborious, and also we could not be sure that the 
values of k chosen would fairly represent all types of be- 
haviour in the circuit. Both these difficulties can be 
avoided by the following procedure : the cubic (14.13) 
has always one negative real root, call this a, and 
suppose that it has also a pair of complex roots db if\ 
of modulus r. Then we know that * the sum of the roots 
of (14.13) equals ~ kb, the sum of their products two at 
a time equals i, and the product of the roots equals k, 
Or, written in full, 



ocr 2 = k, . . . (i4 l6 ) 
where r 2 = 2 + r) 2 . . (14-1?) 

If we regard a as a parameter, and solve (14.14) to (14.16) 
for , r z and k in terms of it, we find 



' 



That is, a, i^/(r 2 2 ), with the values of r 3 
and given by (14.18), are the roots of (14.13) with the 

* Charles Smith, Treatise on Algebra. If a u . ..., oc n are the 
roots of 

x n + a^- 1 + . . . + _!* + a n = o 

then Zoc! = a i} Zoc^ = a z , Sa^ocs = a a> . . ., 

ajoc, . . . a n = (-)*. 
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value of k given by (14-18). If the value of given by 
(14.18) turns out to be greater than r z , the same argument 
shows that there are three real roots, 

- <x, and | -j- ^/ (2 _ r ) 
for the corresponding value of k. 




Naif 5um of roots 



FIG. ii. 

In Fig. 1 1 the position of the complex root $ -\- ii) 
of (14.13) is plotted as a function of the parameter a, 
for the values 2, 4, 8, 9, 10 of &. When <x = o this root is 
i for all values of b, and as <x increases the curves shown 
are traced out. The positions of the roots for the values 
0-2, 0-4, 0-6, i-o, 2-0 of a are marked by successive crosses 
on the curves, k is given by (14.18), and, when there is 

4 
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a complex root, increases as a increases, that is, as we 
move round the curves. If b < 9 there is always a pair 
of complex roots of (14.13). If b 9 there is a real triple 
root for k == i/(3V3) at a = i/Vs ; for all other values 
of k there is a pair of complex roots. If b > 9 there are 
three real roots for some values of k ; these are indicated 
by plotting (dotted in Fig. n) half the sum of the real 
roots other than a as abscissa, and half their difference 
as ordinate, measured below the axis. 
. T^ us a diagram su ch as Fig. n gives an immediate 
indication of the way in which the roots of (14.13) "vary 




FIG. 12. 

with k and b. To study the way in which the current I 
derived from (14.12) varies with k and b, we select typical 
values of b from Fig. 11, for example, the transition value 
b = 9, a value of b less than 9, and a value greater than 9. 
On each of these curves we choose a few representative 
points and calculate I for these using (14.6). 

In Fig. 12 the resulting set of curves for the case b g 
is shown. Here RI is plotted against n^t for the values 
0-05, o% 0-2, i/V3, 2, 4, 8, 16, 32 of a, to which 
correspond the values 0-049, '93> 'i52, 0-192, 0-270, 
0-469 0-901, 1-784, 3-558 of k. Clearly these curves, 
together with two or three such sets for other values oi 

V^J 1 give a fairly com P lete indication of the way in 
which the transient behaviour of the circuit depends on 
k and b. 
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15. Circuits which include vacuum tubes 

Transient problems in valve circuits are easily treated 
by the Laplace transformation method, provided, of course, 
that the valves are operated in linear parts of their char- 
acteristics. Only triodes are considered here, but similar 
considerations apply to other valves. 

The characteristic equation of the triode supplies a 
linear relation between the anode current, I a , and the 
anode, grid, and cathode voltages V a , V,, V , respectively, 
namely p i ffl = (V a - V ) + tfV, - V c ), . (15.1) 
where p is the anode slope resistance, and fj, the amplifica- 
tion factor. In transient problems, the quantities V a , V , 




a 




FIG. 13. 

V c , I appearing in (15.1) are the variable parts of _ the 
voltages and current caused by the transient behaviour 
being studied ; they are superposed on the constant vol- 
tages and anode current caused by the anode, and grid bias, 
batteries. The grid current is negligible in normal cir- 
cumstances, and here is always assumed to be zero. In 
any problem it must be decided whether the interelectrode 
capacities in the valve itself can be neglected if not, 
they are easily included as shown below. By virtue of 
(15.1) it is possible to replace the valve by an equivalent 
circuit element consisting of resistance p and a voltage 
_ ^(V V ) inserted between anode and cathode ; 
this is the common practice, but it is a little more funda- 
mental to work directly from (15.1) as will be done here. 
As a first example we consider the type of problem that 
occurs in the amplification of transients : in the circuit 
of Fig. 13 interelectrode capacities are neglected, and 
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steady conditions prevail due to grid bias and plate batteries 
(not shown in the figure), when at t = o a known addi- 
tional voltage V ff is applied- to the grid ; it is required to 
find the change, l a , in the anode current. Since V ff and 
the changes in currents and voltages it causes are super- 
posed on the steady values, we do not need to know the 
latter, and need only write down equations for the changes, 
all of which vanish at t = o. All the currents and voltages 
marked in Figs. 13-16 are to be understood in this sense. 
We notice that (15.1) implies a convention of sign, namely, 
if V a > V , I is flowing in the direction of the arrow in 
Fig. 13 ; the same current I a also flows from the cathode. 




FIG. 14. 
The subsidiary equations for the impedances z a and 



are 



the signs being determined by the convention stated above. 
Also the transform of (15.1) is 

pi. = V B + /*?. - (x + /OV . (15-4) 
Solving (15.2), (15.3) and (15.4) we find 

{p + * 8 + (i + /*> = /*?,, ta-S) 
and from this l a can be found if V g is known. 

As a more definite problem we take the case of resistance- 
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capacity coupling, Fig. 14, and find the value of V 2 (which 
would be the grid voltage of a second valve) corresponding 
to constant voltage E applied at t o to the grid of the 
first valve. Here we have, in the notation of Fig. 13, 



* 



E 
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FIG. 15. 

where (15.5) has been used in the last line of (15.6). Since 
the denominator of (15.6) is a quadratic in p, it is easy to 
write down V 2 . 

To illustrate the effect of inter-electrode capacity and 
a rather more complicated system, we consider the circuit 
of Fig. 15. Here grid-plate capacity is included in s 3 , 
grid-cathode capacity in z, and anode-cathode capacity 
in * 2 . We seek the anode voltage V a due to voltage V 
applied at t = o at the terminals AB when steady con- 
ditions prevail. 
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Choosing currents as shown in Fig. 15 to give zero grid 
current, the subsidiary equations are 

'*A = ?-? . (15-7) 
*& + 1, - l a ) = V. . . (15-8) 

* 8 (I, - l a ) = V. - V, (15-9) 

.* a l a =-V a , . - (15-10) 

and the transform of (15.1), which is, in this case, 




FIG. 16. 



Solving we find 



V = 

* n. ^^ 



(I5-I2) 



whence V a can be found if the impedances z are known. 
Finally we treat an oscillator circuit (Fig. 16) from the 
present point of view. Here there is a closed system in 
which the grid voltage is not prescribed from outside, but 
is free to adjust itself. As always, there are steady currents 
and voltages due to plate and grid bias batteries, and the 
I,l a , V a , V ff are departures from these : clearly one solution 
is to have all of these zero, and the problem is to find how 
this solution behaves if the system is disturbed. To fix 
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ideas * suppose there is a change of current I in the 
inductance L at t = o, then the subsidiary equations are 

Rl = ? + Li . I . I 



(15.15) 



The equations (15.13) to (15.16) can be solved for 
I, I a) V a , V g , in terms of 1 and the circuit constants. All 
the solutions will have for denominator 



i 



o i 

Cp Cp 

i o 

- i , 



1 , (15.17) 

and the solution contains terms in ef***, where a r , r=i, 2, 3, 
are the zeros of (15.17). If all the oc r have negative real 
parts the disturbance dies out gradually and the system is 
stable ; if any of the oc r has a positive real part the system 
is theoretically unstable, actually it executes an oscillation 
of amplitude determined by the limits to tshe linearity of 
the characteristic. The condition for this, which is the 
condition for the circuit to act as an oscillation generator, 
is obtained by applying Hurwitz's criterion, (12.3), to 
(15.17); it is 
C[p(L + LO + (i + ^(RLt + LRa)] 



1 6. Filter circuits 

The essential feature of such circuits is that a number of 
precisely similar circuit elements are arranged to form a 

* This is done solely as an illustration, all that is needed is 
the determinant A which does not involve the initial disturbance. 
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repeating pattern. Mathematically this implies that we 
can write down a general relation, called a difference 
equation, connecting the transforms of the currents in 
successive elements, and by solving this we get a general 
expression for the current at any point. 

As a definite example which includes many cases of 
practical interest we consider the circuit of Fig. 17. Here 
there are m impedances z, and m i impedances %', 
arranged as shown. Yoltage V is applied at the beginning 
through impedance ss {> and the whole is terminated by 
impedance # . We require the current I r at any point due 
to the voltage V applied at t = o with zero initial conditions.* 



I t _f~* 

\ vi; 



VUf 




T T 

FIG. 17. - m " " 

The subsidiary equations for the successive meshes of 
Fig. 17 are 



The equations (16.2) which connect the transforms of 
any three consecutive currents are a system of difference 
equations. We seek a solution of them of type Ae^, 
where A is independent of r. Substituting this in (16.2) 
gives 

z (2$ -f #')0 e + #e 2& = o, 

/ 
that is cosh Q= i + -. . . (16.4) 

* If the initial conditions are not zero, as in the discharge of 
a. charged filter circuit, there -will be terms on the right hand 
sides of the subsidiary equations which depend on the particular 
impedances z and z' in the circuit. 
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. Thus, if 6 is defined by (16.4), a solution of (16.2), and 
in fact the general solution * O f it, is 

It follows from (16.4) that sinh 2 \Q = #74*, and, since 
both signs of 6 occur in (16.5), the square root with the 
positive sign may be taken, -which gives 

sinh 40 (l\*. . . (16.6) 




-* 



2 2 n 2z 



FIG. 18 (a), FIG. 18 (6). 

To find the constants A and B in (16.5), this must be 
substituted in (16.1) and (16.3), which give two equations 
for A and B, namely 

(*, + * *e)A + (*< + *- %<T*yB = "^ 
A(# -f * - *<r e )e + B(* + a - ar^X"* == o. 
Solving these and substituting in (16.5) gives finally 

= 

sinh (?n r)fl g sinh (OT r 1)0 _ 



) sinh wi^^Csfa+^-f-a^) sinh (wi i)ff+2 a sinh (/w a)0 

(16.7) 

This general expression simplifies greatly when_ ^ 
and % i have values corresponding to the common termina- 
tions. The network of Fig. 17 may be regarded as built up 
of either T elements (Fig. i&z) or ff elements (Fig. 180) 

*The theory of difference equations is analogous to that of 
differential equations. (16.3) corresponds to the c m Pl^ e ta jT 
function, and, if there had been functions of r on the right hand 
side of (16.3), a particular integral would have had to be includetj, 
A brief reference to the theory is given in Piaggio, Differential 
Equations, p. 315. The Laplace transformation method can ^ be 
applied to difference equations, cf. Gardner and Barnes, loc. ctt. 
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connected in. tandem, and the general expression (16.7) 
takes specially simple forms when circuits built up in this 
way are open or short circuited. For example : 

(i) #o = z { o : Fig. 17 with no series impedance, 
short circuited at the far end, 

T V cosh (m-r^ )0 

r zz sinh (m - 1)6 sinh 0' ' ^ Ol0 ' 
(ii) % = oo, ss t o. The same, but open circuit at 
the far end, r?- i. / \/i 

I = V sinh (m r}d 

r ~~ 2z cosh (m - i)0 sinh 0' ' ( 9 ' 
(iii) #o = z f = *'. A circuit of 01, T sections, with 
.no additional series impedance and short circuited at the 
far end. Here, using (16.4) we get 

T V cosh (m r)9 
Lr ~z sinh 6 sinh mB' ' ' (l6 ' IO) 
(iv) * 4 = Jar', # = oo. The same as (iii), but with 
open circuit at the far end, 

j _ "V sinh (m r)9 
r ~ z sinh 8 cosh mff ' ' ( ^ 

( v ) %i = %'> ^ = *' + zx> A circuit of m + i, ff 
sections, with no additional series impedance, and open 
circuited at the far end, 

T = V sinh (-r + i)g / , N 

r 2^ cosh (m + i)0 sinh & ' ( ' I2) 

If "9", z, z', z , z t are known functions of p, (16.4) gives 
cosh 9 as a function of p, and using the formulas 

cosh = a"' 1 f cosh 5 - cos \ fcosh - cos 

\ an/ \ an 



. . . (cosh - cos (2n ~ l)?rN ) (J6.I7) 
\ an / v *' 

sinh w0 = 2"- 1 sinh (cosh cos -\ ("cosh - cos \ 
\ nj\ n J 

. . . (cosh 0- cos ^^, (16.14) 
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it appears that (16.7) can always be expressed as a quotient 
of polynomials in p. Both (16.13) and (16.14) no ^ ^ or 
any positive integer n, even or odd. It is often advisable 
to use them to see the exact form of the denominators of 
the transforms of the currents, though actual evaluations 
of the currents I r are best carried out by working in terms 
of 6 as in the examples below. 

EXAMPLE i. Constant voltage E applied to a circuit 
of m, T sections in which %' R, % i/Cp, with no 
additional series impedance, and open circuit at the far 
end. 

Here (16.4) gives 

cosh 6 = i + |RQ>. . . (16.15) 

And, using (16.11), 

CE sinh (m - r)9 ^ 



^ 
T sinh cosh mQ 

To evaluate I r we need to know the zeros of the de- 
nominator of (16.16). These can be found from (16.13) 
with the value (16.15) of cosh 6, or by noticing that cosh mQ 
is zero when 






~ Taking more values of s than those in (16.17) would 
result in counting the same zeros over again. By (16.14), 
sinh 9 is a common factor of the numerator and de- 
nominator of (16.16) so 9 = o need not be considered 
in using Theorem II. 

The values of p corresponding to the values of 6 given 
by (16.17) are > using (16.15), 

5 = , i, . . ., OT-I. (16.18) 



To evaluate I r it is best to use the form (2.8) of Theorem 
II, and for this we need 
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d 



r d 

I j- (sinh cosh mff) 

= -^ (sinh cosh m0) . -5- 

\_dd v dpjp^a,, 

cosh cosh md -\- m sinh sinh 



L* 



sinh 



sinh 
i)tn- 



) ...... (16.19) 

Then, applying (2.8) to (16.16), using (16.19), & 

, sE^=J . ( OT - r )(2s 

l r = - > ( V e*s* sin ^ - ^ - ! 






EXAMPLE 2. Voltage V = sin tat is applied at t = o 
to a filter of m, T sections, in which *' = L#>, x = i/Cp, 
with open circuit at the far end. 

Here (16.4) gives 

cosh = i + LQ> 2 .. . . (16.21) 
And from (16.11) 

T _ Cpa) sinh (m y)0 
r ~(> a + co a ) sinh cosh mff ' ( l6 ' s2 ) 

The zeros of the denominator of (16.22) are + ico, 
and the values of p corresponding to 
a _ ( 2S + i)*w 

^ - ' ' = ' *> ' ' m ~ * (16.23) 

Using (16.21) these are seen to be 

P = ip,, . .... (16.34) 
where = \* f i - c &Ltl 

r cos ^r 



and it is assumed that none of these equals to. 
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Proceeding as in Example i we find 

-r- (sinh 6 cosh 0)T = ( ) a+1 LCwz/3 3 . (16.26) 

J*P Jp^ip* 

And thus from (2.8) and (16.22) 

I r = -=- "S* s ==- cos fl.i cos ^ 

mL *~ o> 2 B z a rs -2m 

t~ o / y 

, C Ceo sinh (m -rWzty) V1 . . _, . ,. "I / x % 
+ 1 . , fl/ . \ , L- \ g * -f- Conjugate V (16.27) 
U sinh 0(to>) cosh i0(ia>) J te J 

where 0(zo>) is the value of 6 corresponding to p = uo. 

The second term of (16.27) * s tne steady state part of 
period 27r/co. To evaluate it, we see from (16.21) that 

cosh 0(w) = i JLCoj 2 , . (16.28) 

and this gives rise to two possibilities : 

(i) if o> 2 < 4/LC, 0(to>) = z'S, 

where 8 = cos- 1 (i -|LCa> 2 ), 

and in this case the part of L of period 27r/co is 

~ sin (wz r) S ^ , ., , x 

Ceo . \, ~ cos o>t. . . (16.29) 

sin 8 cos mo 

(ii) if a; 2 > 4/LC, 0(co) = ZTT -(- A, 

where A = cosh- 1 (LCa> 2 i), 

and in this case the part of I r of period 277/co is 

, N 1r , sinh (m r)A / u \ 

( ) r - 1 Coj . . : = ~ cos cot . . (16.30) 

v ' sinh A cosh wzA \ J J 

As r is increased, (16.30) decreases exponentially, while 
on the other hand (16.29) oscillates. That is, if to 2 > 4/LC 
the term of period 27r/co dies out exponentially as we move 
along the circuit, while if to 2 < 4/LC this term oscillates 
trigonometrically. The circuit is, of course, a low pass 
filter, and we have found both its transient and steady state 
behaviour for applied voltage sin cof. 

17. Vibrations of mechanical systems 

Problems on the motion of bodies of concentrated inertia, 
connected by shafts or springs of negligible inertia, lead 
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very readily to systems of ordinary differential equations 
which have to be solved with given initial conditions. In 
order tp make these differential equations linear, the dis- 
sipative forces are restricted to the case of resistance to 
motion proportional to the velocity, though in some 
problems Coulomb friction (constant and opposite to the 
direction of motion) can be treated. Rotating systems 
only will be considered here, since 'they are of perhaps 
a little more interest than the corresponding problems in 
linear motion, and lead to precisely similar equations. 

We write throughout this section, J for the moment of 
inertia of a wheel about its axis, 9 for its angular displace- 
ment from a standard position, D for differentiation with 
respect to time, ^ = D9 _ > .^ (l?>l) 

for angular velocity, kca for resistance to motion proportional 
to angular velocity, T for torque, and A for the stiffness of 
a shaft, that is the torque required to produce unit relative 
angular displacement between its ends. 

The equation of motion of a wheel of moment of inertia 
J, acted on by torque T, and with resistance to motion 

ka) ' is (JD-M)o, = T. . . (17.2) 

Also the torque T required to produce a relative dis- 
placement 6 between the ends of a shaft of stiffness A is 

T = A0. . .. . (17.3) 

Using these equations we can write down differential 
equations for the motion of any combination of wheels 
connected by shafts of known stiffness. A number of 
examples of the differential equations which arise in such 
problems is given below : clearly these can be solved just 
as in 5 with given initial values of the o> and 9. 

EXAMPLE i. A shaft of stiffness A has one end fixed, 
and on the other a wheel of moment of inertia J (Fig. 19(0)). 
Torque T, a given function of time, is applied to the wheel, 
and there is resistance ku> to its motion. 

If 9 is the displacement of the wheel from its equilibrium 
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position, the retarding torque due to twist of the shaft is A0, 
and the differential equation for the motion of the wheel is 

(JD + *) = T - A0. . . (174) 



ir 



T 

- 



FIG. 19 (a). 



FIG. 19 (&). 



Equations (17.1) and (17.4) can be solved with given 

initial values of 6 and co in the usual way. If preferred, 

CD can be eliminated ; this gives the second order equation 

(JD 2 + 6D + A)0 - T. . . (17-5) 
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FIG. 20 (6). 



EXAMPLE 2. Three wheels are connected by shafts of 
stiffness J\ and A x as shown in Fig. 20(0), in which the 
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symbols used all have the same meanings as before, 6, # a , a 
being displacements of the wheels from a reference position 
in which the shafts are unstrained. Torque T is applied 
to the wheel J. 

Here the differential equations for the motions of the 
wheels J, J 15 J 2 , respectively, are 



CJiD + ^K = - Aft - 6) + M*. 
(J^ + ^K^-A^-^), . 
and a) = D0, ^ = D^, a> z = D0 2 - . 



(17.7) 

(17.8) 

(17.9) 



This system can be solved in the usual way with given 
initial values of 6, 8^ 2> u>, o^, o 2 . Clearly the problem 
of any number of wheels connected by shafts of given 
stiffnesses can be solved in the same way ; if all the wheels 
have the same moment of inertia and the same resistance 
to motion, and all the shafts the same stiffness, the sub- 
sidiary equations form a chain of the type met hi 16 
and are treated in the same way. 

EXAMPLE 3. ^4 geared system. 

First it is necessary to study the type of complication 
introduced into the system by 
gearing. Suppose a gear of 
radius r drives a gear of radius 
r a> tne i r angular velocities being 
wj and cu 2 in the directions of 
the arrows in Fig. 21. Then 
if #! and a are the angular dis- 
placements of the two gears 
from their original positions 




p* 



(17-") 

At the point of contact of 
the wheels there is a driving 
FIG. 21. force P on the wheel r a , and 
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an equal, oppositely directed, reaction force P on the 
wheel r v Then if T 2 is the driving torque on the wheel r z , 
and T! is the reaction torque on the wheel ^ 



and thus 



(17.12) 




FIG. 32 




FIG. 22 (6). 

(17.11) and (17.12) form a pair of linear connections 
between the angular velocities and torques in the two 
shafts. More complicated systems, such as differential 
gears, give rise to linear connections in the same way. 

As an example we consider the system of Fig. 22 (a) 
in which torque T is applied to the wheel J which drives 
a gear of radius r x and moment of inertia J^ through a 

5 
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shaft of stiffness A. The gear Jj meshes with a gear of 
radius r z and moment of inertia J 2 , which drives a wheel 
of moment of inertia J s through a shaft of stiffness A r 
With T! and T 2 defined as in (17.12) the equations of 
motion are 

T + A^-fl) . (17.13) 
--A(0 l -.0)--T 1 . (17.14) 
(J 2 D + > 2 = Ai(0 8 - a ) + T 2 . (17.15) 
(J 3 D + > 3 = -A^ - 2 ). . (17.16) 

These equations, together with (17.1), (17.10), (17.11), 
(17.12) are sufficient to give the motion for any initial" 
conditions consistent * with (17.10) and (17.11). 

In many cases it is possible to set up an electric circuit 
problem which corresponds f to a given mechanical 
vibration problem, in the sense that both lead to precisely 
the same systems of differential equations. This does not 
make the problem any easier to solve, but it gives an inter- 
esting physical insight into it, also the corresponding 
electrical system may be set up and studied in the laboratory, 
and finally, since the theory of electric circuits is more 
highly developed than that of mechanical systems, it may 
be useful in suggesting further developments in that 
field. 

To see the sort of correspondence J which exists between 
mechanical and circuit problems, compare equations 
(17.2), (17.3), (17.1), respectively, from which the differ- 
ential equations for mechanical problems were built up, 
with the fundamental equations from which the differ- 

* Initial conditions inconsistent with these correspond to 
forcing into mesh at t o a pair of gears inconsistently moving 
or displaced. 

t For a large number of such equivalent problems, cf. Olson, 
Dynamical analogies (van Nostraoad, 1943). 

J The correspondence described here is not the only one ; 
thus, instead of (17.17) we might have used the equation for 
a leaky condenser, namely (CD + G)V = I, and in this way 
another, equally useful, system arises ; cf. Miles, Journ. Acoust. 
Sac. America, 14 (1943), 183. 
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ential equations of circuit problems are built up, namely : 
(LD-|-R)I = V, . . (17.17) 
for an inductive resistance with applied voltage, V ; 



for a condenser with applied voltage, V ; and 

DQ = I, . . (17.19) 

for the connection between charge and current. It appears 
that (17.2), (17.3), (17.1). a nd (17.17) to (17.19) are the 
same sets of equations in terms of different quantities, 
the correspondence between these being set out in the 
Table below, in which the corresponding quantities for 
linear motion of a mass are added for completeness. 

Electrical. Rotational. Linear Motion. 

Charge, Q. Displacement, 6. Displacement. 

Current, I. Angular velocity, w. Linear velocity. 

Inductance, L. Moment of inertia, J. Mass. 

Resistance, R. Damping coefficient, k. Damping coefficient. 

Capacity, C. Reciprocal stiffness, i /A. Reciprocal stiffness. 

Voltage, V. Torque, T. Force. 

Using this correspondence we can now set up electric 
circuits which have the same differential equations as 
the mechanical examples considered earlier. For example 
the circuit of Fig. 19 (b) with inductance J, resistance k, 
and capacity i/A, with current o>, charge 9 on the con- 
denser, and applied voltage T, gives the same differential 
equation (17.4) as the mechanical problem of Fig. 19 (a). 

Similarly Fig. 20 (b) is the equivalent circuit * to Fig. 
20 (a), so that the differential equations found by applying 
Kirchhoff's second law to the closed circuits AA 1 B 1 B, 
AjAaBjjB!, and A 2 A 3 B 3 B 2 , are (17.6) to (17.8) respectively, 
where 6, 6^ 6 Z are the amounts of charge which have 
passed through the inductances J, J 1; J 2 at time t, so that 

* Leaky condensers in Fig. zo(b) correspond to resistance to 
motion proportional to the relative velocities of the wheels in 
Fig. 20 (a). 
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0J and X 2 are the charges on the condensers 
i /A and i/A x respectively. 

Finally we consider the equivalent of the geared system 
Fig. 22 (a). The linear connections (17.11) and (17.12), 

<DI To r 
i.e. _i = _2 = -?, . . (17.20) 

COo, " 1 ' i 

become in the electrical case 

r = ^=7> (*?) 

^2 v l M. 

which are the connections between currents and voltages 
in an ideal * transformer of r 2 secondary and r x primary 
turns. On this understanding the equivalent circuit to 
Fig. 22 (a} is Fig. 22 (b), and it is easy to check that these 
lead to the same differential equations. 

1 8. Servomechanisms 

The ideas involved f may be illustrated by considering 
a rotation control system in wmch a massive body of 
moment of inertia J has to be driven so that its displace- 
ment follows as closely as possible the prescribed dis- 
placement 9i(t) of a pointer. To do this, voltage f 
is fed into an amplifier-motor system which applies torque 
k(6t ) to the body, where k is a constant. This torque 
tends always to bring the " output displacement " 
into agreement with the " input displacement " i5 and is 
proportional to the " error " = f . To provide 
damping for the system, it is usual to " feed-back " into 
the amplifier a voltage depending on the output velocity 
D0 , or on the error velocity D0. The simplest arrange- 
ment is to mount a generator on the output shaft which 
provides voltage k^Dd^ which is fed back to the amplifier 

* Cf. Starr, Electric Circuits and Wave Filters (Pitman, 1934), 
p. 128. 

t Cf. Hazen, " Theory of servo-mechanisms ", Journ. Franklin 
Inst., 218 (1934), 279 ; Minorsky, " Control Problems ", ibid., 
232 (1941), 451, 519. 
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with negative sign. The torque supplied by the motor is 
then k(6 t 6 - &iD0 ), and the equation of motion of 
the body is JD 2 = k(e t - 8 Q - &,D0 ). . (18.1) 

If the motion is started at t = o from rest at zero dis- 
placement, the subsidiary equation is 

kfa^kei, . . (18.2) 






where co 2 = fc/J. ... . (18.4) 

Suppose the input is given a uniform motion 9 t = Q, t , 
starting at t o, so that 



then (18.3) gives 



O feff , 0(^4- 
~ 2 ~"~ 2 + 
It follows that 



cos ni + -(^ 2 ^> 2 - i) sin nt\ (18.7) 



where w 2 = oj 2 l^ 2 ^ 4 . ( 

It is seen from (18.7) that the error 6 9 6 t in this 
case is composed of a constant lag, k]Q , together with 
transient terms. 

If instead of feedback ,p0 , the feedback is arranged 
to be proportional to the error velocity, kjDd, the equation 
of motion is JD 2 = k(- d - k^Dd) . . (18.8) 
or JD 2 + kk^pe + kd = - JD 2 ^. 

If the motion is started from rest at zero displacement, 
the subsidiary equation is _ 

(p 2 + k^p + co 2 )0"= - p z e t , . (18.9) 
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where w 2 = &/J. If, as before, 8 t Q,t, this gives 
x -O 



and 0= --*-**!"* sin *. . . (18.10) 

w N 

In this case the error (18.10) consists of transient terms 
only and there is no permanent lag as in (18.7). 

In practice the voltage kjl)6 or kj^d is usually not 
fed back directly to the, amplifier, but applied to some 
circuit, and voltage is tapped off from some part of this 
circuit and fed to the amplifier. In this way rather more 
complicated transforms arise. 

Another complication of practical importance is the 
presence of time lag in some part of the system. To 
illustrate the effect of this suppose that, in the first case 
considered above, the voltage 6 t 9 Q k^Dd^ instead of 
being fed directly to the amplifier, is applied to an induc- 
tance L and resistance R in series, and that the voltage 
across the resistance is fed to the amplifier as before. 
Here there is a time delay determined by the time constant 
of the L, R system. The subsidiary equation for the 
current I in the L, R circuit, with zero initial conditions is 

(Lp + R)I = 6 t - (i + krffa, . (i8.ii) 
and the voltage V across the resistance R is given by 

V = RI. . . . (i8.ia) 

Also the subsidiary equation for the motion of the body 
of moment of inertia J, with applied torque KV is 



= kV. . . . (18.13) 
From (i8.ii) to (18.13) it: follows that 



where o> 2 = kj], c^ = R/L. . . (18.15) 

The error may be evaluated, for given <) . when the 
roots of the period equation 



7* 

are known,. Hurwitz's criterion (12.3) for the roots of 
this equation to have negative real parts, that is for the 
system to be stable, requires 

^>iK. . . . (18.17) 

Thus if there is time delay in the system it is necessary 
for the amount of feedback to be greater than the minimum 
amount (18.17) or tn ^ e system will not be stable. The 
way in which the roots of (18.16) and the values of 
vary as k^ is increased may conveniently be studied by the 
method of 14. 

19. The impulsive function 

In dynamics the idea of a very large force acting for 
a very short time and such that the time integral of the 
force is finite is extremely useful, and the "impulsive 
motion " of dynamical systems is extensively studied. 
In circuit theory the ideas of an impulsive voltage (that is, 
a very great voltage applied for a very short time and such 
that its time integral is finite) or an impulsive current, 
are equally fruitful.* 

It is convenient to have a symbolic representation of 
such a function and this may be obtained as follows. 

Consider the function </>() defined f by 




(19.1) 
= o t >e 

* For an interesting account of impulsive functions and their 
applications see van der Pol, jfourn. Inst. Elect. Eng., 81 (1937), 
381. 

t Many other functions have the same properties, e.g. the 
triangular function, Fig. 23 (&), or the function 



Fig. 23 (c), except that the latter has its maximum at t ~ o. 
It is more usual to choose a function symmetrical about t = o, 
such as Fig. 23 (c), but there is some advantage in the present 
applications in being able to think of, say, an applied voltage 
which is zero at t = o and becomes very large immediately after 
this instant. 
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where e is very small. This is very large in a very small 
region to the right of the origin (Fig. 23 (#)) and zero 





f a+e 



O 2 

FIG. 23 (a). FIG. 23 (5). FIG. 23 (c), 

elsewhere, and, no matter how small e is, 

f , \ 

<j)(f)dt = i. . . . (19-2) 
J-oo 

Also, if f(t) is a continuous function 

TOO j /-<H 

J-oo eja 

= /(*) + O(e), - (I9-3) 

where O(e) is written for a term involving e as a factor, 
and thus which tends to zero as e -^- o. 

We now define the " impulsive " function 8(i) to possess 
the properties of (j>(t) as e ^ o that is : 

S(t) = o for t < o, and for t > o, except for a"| . v 
vanishingly small region to the right of t = o /' 

poo 

8(t)dt =i, . . . (1.9.5) 

J-OO 

fco 

f(t}^(t a)dt=-f(a\ . . (19.6) 
J oo 

if /(i) is a continuous function. 
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It follows from (19.6) that 

i, . . . (19.7) 



o 

that is, the Laplace transform of the 8 function is unity. 
We may now regard a blow of impulse P as a force 
PS() ; an impulsive voltage whose time integral is E, 
as a voltage ES(i) ; and an impulsive current which transfers 
charge Q as a current QS(), and proceed to solve problems 
involving them in the usual way : the validity of the 
solutions so obtained is discussed later. 

EXAMPLE i. A mass M attached to a spring of stiffness 
oj 2 M is set in motion at t = o from rest in its equilibrium 
position by a blow of impulse P. The equation of 
motion is 



and, using (19.7), the subsidiary equation is 

M( 2 + o) a ) = P. 
Thus x = (P/Mco) sin tat. . . (19.8) 

EXAMPLE 2. An impulsive voltage ES(i) is applied at 
t = o to an L, R, C circuit with zero initial conditions. 
Here the subsidiary equation (8.6) is, using (19.7), and the 
notation (9.3), * _ Ej> 



Thus, in the case n 2 > o, 

E 

I = _ e -A*<{ w cos nt ju, sin nt}. . (19.9) 
Ljfi 

The solutions of Examples i and z above have been 
obtained by the ordinary Laplace transformation procedure, 
using (19.7) for the transform of S(), and, since any use 
of the 8 function must be regarded as symbolical only,* 

* The S function was denned as the limit of some function as 
e -> o. Any mathematical operations on this function, such as 
differentiation or integration, would involve the interchange of 
order of two limiting processes and be very difficult to justify. 
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we must inquire whether these solutions are mathematically 
justifiable. The simplest way to do this is to verify that 
they do in fact satisfy the differential equations and initial 
conditions of their problems. Consider the solution 
(19.8) from this point of view: it satisfies the differential 
equation MD 2 # + Mco 2 * = o, *> o, . (19.10) 



for the motion of the mass M for t > o ; also as t -+ o 
in (19.8), x -> o, and Dx -> P/M. Now P/M is the velocity 
with which, according to the theory of impulsive motion 
in dynamics, the mass M would be started by a blow of 
impulse P. Thus the solution obtained by the formal 
Laplace transformation procedure applied to Example i 
does in fact giye the complete solution of the dynamical 
problem correctly, and in addition it avoids the necessity 
of evaluating the initial velocity due to the blow. It can 
easily be verified * that the same is true for a general 
dynamical system of n degrees of freedom set in motion 
by blows. 

Precisely the same is true for the electrical case in 
Example 2. Here, as t -> o the value of I given by 
(19.9) tends to E/L, and we need an electrical analogue 
of the theory of impulsive motion in dynamics to show 
that an impulsive voltage ES(), applied to an L, R, C 
circuit in which no current is flowing, will instantaneously 
produce a current of this magnitude. This can be deduced 
from the differential equation for current and charge in 
the circuit, namely 

(LD + R)I + = E8(). (I9-") 

Proceeding as in dynamical theory, we integrate this 
over a very small time, T : the integral of S(i) over this 
time is i ; that of the finite quantities I and Q is very 
small ; and that of DI is the difference between the final 

* By an easy extension of the method in Carslaw and Jaeger, 
loc. cit., 35. Jaeger, Phil, Mag. (7), 36 (1945), 644, gives a 
complete verification for all the types of problem arising in this 
section and 20. 
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value, ly., and the initial value (zero) of the current I. 
Thus LI, = E, . . . (19.12) 

and the theory of impulsive motion gives the current 
I/ = E/L, to which the solution obtained by the Laplace 
transformation procedure was found to tend as t -> o. 
It is easy to verify that (19.9) satisfies the differential 
equation of the problem for t > o, so it does in fact give 
the correct solution of the physical problem. 

It can be proved in general that for a complicated circuit 
with various impulsive applied voltages, the solutions 
obtained by the Laplace transformation procedure, treating 
B(t) as having (19.7) for transform, do satisfy the equations 
of circuit theory, and, as t > o, the currents tend to the 
values which would be found from the equations of 
impulsive motion. 

20. Justification of the solutions of electrical network 
problems 

It was remarked in 7 that, for a system of linear 
differential equations, a general verification could be given 
that the solutions obtained by the Laplace transformation 
method satisfied the differential equations and initial 
conditions, provided a certain determinant did not vanish. 
The problems of circuit theory often give rise to a system 
of differential equations to which this verification applies 
immediately. In general, however, they give a combina- 
tion of differential or algebraic equations (10.1) and alge- 
braic equations (10.6) ; it may happen that the specified 
initial values of the currents and charges do not satisfy 
some of these equations, and in that case there will be an 
impulsive readjustment of these initial values. It is in 
fact true * in all cases that the solutions obtained by the 
Laplace transformation method, used as in 10, giving 
each inductance its own initial current, and each con- 
denser its own initial charge, in the subsidiary equations, 
do in fact give solutions of the equations of the problems 

* Jaeger, loc. cit. 
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for t > o, which as t -$ o tend to the values which exist 
in the circuit after the impulsive readjustment of initial 
conditions, if this takes place. The general theory will 
not be given here, but only examples illustrating the 
different cases which can arise. 

EXAMPLE i. Steady current E/R is flowing in the 
circuit of Fig. 24 with the switch S closed. At t = o 
the switch is opened. 



innf^- 
L, i, 



FIG. 24. 

Since the initial current is E/R in L, and zero in 
and L 2 , the subsidiary equations are 




. (20.1) 



Solving we find 



Rp R{(LL! + LL 2 + 
And thus 



+ 



H- L.,)} 



(20-2) 



As 



E(LL X -|- LL 2 ) 
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and, since this does not equal E/R, it appears that there 
has been an impulsive readjustment of current between 
the three inductances just after the switch was opened. 

To investigate what happens in this case, and to verify 
this result, consider the equations of the system for t > o, 
namely 

(LD + R)I 4- L 2 DI a - E, . . (20.5) 

L 2 DI 2 - LjDIi = o, . . (20.6) 

I a + l!=I. . . (20.7) 

The values E/R, o, o, which the currents I, I 1} I 2 had 
before the switch S was opened, do not satisfy (20.7). 
Suppose that these change to I<), 1^) and I 2 () in a very 
short time r following the opening of S at } = o, then, 
integrating (20.5) and (20.6) from t = o to =r,'and 
neglecting the integrals of the finite quantities E, I, l lt I 2 
over this very small time, we find 

L{I<> - (E/R)} + L 2 I 2 C> = o, . (20.8) 

L 2 I 2 W - L^W = o. . (20.9) 
Also I<), I^o) and I 2 () satisfy (20.7), that is 

KO) = ^(0) + yo). m . ( 20 . 10 ) 

Solving (20.8) to (20.10) gives the new values of the 
currents, and, in particular, I<) is found to have the value 
(20.4). If the voltage drops across the inductances were 
evaluated in the usual way, they would be found to contain 
S functions, that is, the finite changes of current may be 
regarded as being produced by impulsive voltages. 

EXAMPLE 2. The condenser C in the circuit of Fig. 25 
(p. 78) is charged to potential E, and the condenser Q is un- 
4 charged, when at t o the switch S is closed. 
Here the subsidiary equations are 



(20.11) 
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Solving we find 

. _ 



R(C -f- 
_ _ ECQ ' EC 2 



Thus I=- 



ECCj 



R(C 



- (30.12 



C ^ 


v v\r~ 




11 ^ ' 


I- 1| 


*! 


He " ' \ 


r 1 ' 
11 




< "c, 



FIG. 25. 

The term [ ECC^C -f C^)]8(t) in (20.12) represents 
an impulsive current which instantaneously exchanges 
charge ECC 1 /(C -(- Qi) between the two condensers ; 
the second term represents the subsequent discharge of 
the condensers through the resistance. Whenever a 
network includes a closed circuit containing condensers 
only, there will be an instantaneous exchange of charge in, 
this way between the condensers if these are incompatibly 
charged. If we had used the subsidiary equations for 
charge instead of those for current, this finite change of 
charge would have been found (in the same way as the > 
finite change of current in Example i) but as we have 
worked with current, the impulsive current which transfers 
this amount of charge has appeared. The result (20.12) 
can be verified by writing down the differential equations 
for the charges on the condensers and proceeding as in 
Example i. 
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EXAMPLE 3. Steady current E/^ is flowing in the 
primary of a transformer, Fig. 26, and the secondary 




1 



M 



FIG. 26. 



current is zero, when at it = o the switch S is opened. 
Writing Ij and I 2 for the primary and secondary currents, 
the subsidiary equation for the secondary is 

Mpl x + (Lap + R 2 )I 2 = ME/Rj . (20.13) 
Also L = o, for t > o, 



so 



= o. 



Thus- from (20.13) I a - 






K . . (20.14) 

Here as t -> o, I 2 -> ME/RjL^ There was a sudden 
change in the currents when the switch was opened 
which could have been studied as in Example i. 

EXAMPLE 4. The problem of n, Example i, for the 
case of L]L 2 M 2 (perfect coupling). . Here (11.2) becomes 

- ME . ' , 

In : 



and thus 



ME 



. (20.16) 

It is known from general theory that LjL 2 > M 2 ; the 
limiting case L]L 2 = M 2 is of interest for two reasons : 
firstly it may be very v nearly approached in practice ; 
secondly, if L X L 2 = M 2 , the degree of the denominator 
in 1 2 is reduced by one, here a quadratic has become linear 
and similarly a cubic becomes a quadratic, and the algebra 
is greatly simplified. 
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The significance of perfect coupling from the theoretical 
point of view * is best seen by considering the differential 
equations of the circuit of Fig. 26 : these are 

(LJ> + ROIi + MDI 2 = E 

( 7) 



If LjL.2 = M 2 , multiplying the first of (20.17) by L 2 , 
the second by M, and subtracting, we obtain 

L a Rj! - MR a I a = LaE. . (20.18) 

Thus in this case we can form an algebraic equation by 
elimination from the given differential equations, and if 
the given initial currents do not satisfy (20.18) they will 
change suddenly to hew values which can be calculated 
as in Example i. 

The examples above illustrate the types of problem in 
which there may be sudden changes of charges or currents 
on closing switches, and the way in which these could 
have been found has been sketched briefly. A general 
discussion along these lines can be given, which shows 
that the Laplace transformation procedure, applied in the 
usual way, does in fact always give the correct solution of 
the physical problem. 

/ 
EXAMPLES ON CHAPTER II 

1. Alternating voltage E sin (nt -f- e) is applied at t = o to 
inductance L and capacity C in series. n z = i/LC. If the 
initial current and charge are zero, show, using (i.u) and (1,12), 
that the current at any time is 

(E/2ttL){n sin (nt + e) + sin e sin nt}. 

2. Constant current (E/R) is supplied for f > o to a parallel 
combination of inductive resistance L, R and capacity C with 
zero initial conditions. Show that the voltage drop across the 
condenser is 

E - Ear*" {cos nt + [(R/anL) - (i/RC)] sin nt}, 

* This is an example of the exceptional case referred to in 
7, since the determinant discussed there is LjL 2 M 2 for 
the system of equations (20.17). 
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where /* and n are defined in (9,3). Deduce the solution of 
9, Ex. 3. 

3. ,In the problem of n, Ex. 2, show that the charge on the 
first condenser is 



(CE/aVs){(i + VS) cos * 

- (i - VS) cos wi[i(3 + 

4. Regarding an imperfect condenser as capacity C in series 
with resistance R 1} the whole being shunted by conductance G, 
show that the subsidiary equation for voltage V applied to such 
a condenser in series with an inductive resistance L, R is 

_-r , T ? Q 

- V+ LI ~ 



5. In the bridge of Fig. 9, z\ consists of resistance Rj. and 
capacity C x in parallel, z z consists of resistance R s and capacity 
C 3 in parallel, and st a and # 4 consist of resistances R a and R 4 , 
respectively. Show that there is a balance for -all types of applied 
voltage if CI/CB = R4/Ra = RB/RI> 

6. In the circuit of Fig. 6, LaL 2 = M a , and steady-state con- 
ditions prevail due to alternating voltage of peak value E and 
frequency W/ZTT in the primary. At the instant of zero secondary 
current the condenser is instantaneously discharged, show that 
at time t after this instant the secondary current is 



R a C a 

7. In the resistance-capacity coupled amplifier of Figs. 13, 
14, with # = o, show that if the grid voltage is changed by 
E at t o, the change in grid voltage of the next valve is 

V a *= - 

8. A mass M hangs by a spring of stiffness M a . At t = o, 
when the mass is at rest in its equilibrium position, the other end 
of the spring is given a downwards motion a sin a>t. Show that 
the displacement of the mass from its equilibrium position at 
time t is nn 

Uft - . -i 

r - ;[w sin nt n sin eat]. 
w a n 4 

9. Constant|voltage E is applied at t = o to a condenser of 
capacity Cj in series with a leaky condenser of capacity C and 
leakage conductance G. Show that, if the initial charges are 
zero, the current is 
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21. Heaviside's shifting theorem 

THEOREM VI. // y(p) is the transform of y(t), and 
a > o, then e~ ap y(p] is the transform of 

y(t a)H(* a), . . (ai.i) 
where * H(t) = o, t < o\ 

= i, t> of 



. (21.2) 




y(t-aV H(t-a) 




FIG. 27 (a). 

Since H(i a) is zero if i < a and unity if t > fl, the 
function y(t a)H(i a) is zero for t < a, and for > a, 
its graph is exactly that of y(t) for f > o, shifted to the 
right through a distance a, as in Fig. 27 (tf) : hence the 
name " shifting theorem ". 

* This is Heaviside's unit function; the notation (21.2) will 
be found very useful in many problems. It has not appeared 
earlier in the present treatment, since problems have been stated 
hitherto in the, form of, say, " unit voltage applied to a circuit 
at t = o ", instead of " voltage H(f) applied to a circuit " : it 
must be noticed that these statements are identical only for 
t > o ; for t < o the first does not specify the voltage but the 
second makes it zero. 

82 
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To prove (21.1) we have only to notice that 

poo foo (co 

erv*y(t a)H(t a)dt=\ er**y(ta)dt=\- e~^e- 

JO Ja Jo 



Using this theorem we can build up the transforms of 
a large number of useful step and broken functions. 



O. 

FIG. 37 x-,- 

EXAMPLE i. y(f) = E, o < t < ff\ 

o, t > a I ' 



= o, t > a J 
i.e. XO = E{i - H(* - a)}, cf. Fig. 27 (6), 

J" = !(* -*-*) 



(21-3) 




7W 

FIG. 27 (c). 

EXAMPLE 2. A half sine wave, Fig. 27 (c), 
y = E sin a>t, o <t < 7r/o>\ 

= 0, * > 7T/W /' 

i.e. y = E sin wt + E sin cj( 

jy 



Eo> 
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EXAMPLE 3. A " square " wave, Fig. 27 

y = E, rT < t < (r + i)T \ 

= - E, (r + i)T < * < (r + a)TJ 

i.e. 3; = E{H() - 2H(* - T) + 2H(i - 2 



o 



y = 



3T 



FIG. 27 (d). 



p 
= - tanh 



T 2T 

FIG. 27 (e). 

EXAMPLE 4. A " saw-tooth " wave, Fig. 
y=P k(t rT), rT < t < (r + i)T, a 
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EXAMPLE 5. A function f(t) of period T, so that 
f(t + rT) = /(), if r is any integer. 

_ foo pT f2T 

f(P) = e~#f(t)dt = e-Ptf(t)dt + e 
Jo Jo JT 



= f *-**/(*>& + <r* T I* e- 
Jo Jo 



(21-5) and (21.6) can of course be deduced from this result. 

EXAMPLE 6. The transforms of the fun ctions in Examples 
3 to 5 have all involved hyperbolic functions. A technique 
which is of great importance both here and in connection 
with partial differential equations, is that for finding the 
function which has a given transform of this type. To 
do this we expand the transform by the Binomial Theorem 
in a series of negative exponentials : for example if 

Ei _ e~ 2 ^ T 



.2V 1 e X 1 e ~r e 



_ E 

~P'' 
_ E 

~^' 

Then, by Theorem VI, 
y(t) = Ei 2E(* 

+ 2 E(* - 4T)H(i - 4T) - . . ., (21.8) 
the function y is shown in Fig. 27 (/) (p. 86). 

EXAMPLE 7. A battery of voltage E is applied at t = o 
to an L, R, C circuit with zero initial conditions, and at 
t = T the battery is short-circuited. 



86 
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Here the voltage in the circuit is E EH(i T), and 
by Example i its transform is 



Then, as in (9.2), the transform of the current is 

- - 

and, by Theorem VI, 



I = -{e-# sin nt - -X*- 

TlLj 

for the case 2 > o. 



T) sin n(t - T)}, (21.9) 



2T 4T 

FIG. 27 (/). 



The solution for applied voltages such as those of 
Examples 3 and 4 is found in the same way (cf. Example 5 
at the end of this Chapter). An alternative form for the 
solution in these cases can be obtained by using the ex- 
pansion formula of 33. 

EXAMPLE 8. To find the deflection y at any point of a 
light uniform beam of length /, which is clamped hori- 
zontally at the same level at its ends, x = o and x = l, 
and carries a concentrated load W at x = a. 

Using the S-function defined in 19, a concentrated 
load W at x = a may be treated as a distributed load 
WS(# a). Then, proceeding as in 6, we have to solve 
the differential equation 



EID*y = W8( - a), 

with y = Dy = o, when x = o and x = I. 



(21,10) 
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The subsidiary equation is, by (19.7) and Theorem VI, 



and thus, from (21.11) and Theorem VI, 



- 6EI - 

where 3> 2 an l ^s are unknown constants to be found from. 
the conditions y = Dy =* o, at * = Z. Solving for these 
we get finally, 



6/ 3 EI 

22. Heaviside's series expansion 

This is a useful method for determining the behaviour 
of the solution for small values of the time the solutions 
obtained previously are equally suitable for all values or 
the time : those of this section are useful only for small 
values of the time, but are much more quickly obtained 
than the general solutions. 

Most of the transforms in Chapters I and II have been 
found to take the form 



where g(p) is a polynomial in p of degree n, and f(J>) a 
polynomial of lower degree. If we divide numerator and 
denominator of (22.!) by p n , and expand in ascending 
powers of ifp by the Binomial Theorem, we find an ex- 
pression of the form 

. 5 + ++ . . . (22.2) 



( 

and thus y = % + oj + a^ + . . . . \' 2 -' 2 "^) 

by (i-S), and y is obtained as a series of ascending powers 
oft. 
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EXAMPLE i. In 9, Example i, we had 



i/c 



R \ \ 

TV ; 

EXAMPLE 2. In the problem of 14, (14.1) gives 



Therefore, I = - 

If we letp -> oo in, (32.2), we find 



and since, by (22-3), ! is the value of y when -> o, we 
have proved, /or the form (22.1) o/y, the following Theorem. 

THEOREM VII. Ify(p) is the transform ofy(t), then 

lim y(t) = lim py(p). . . . (22.4) 

t-*0 p-^-ao 

This result is in fact true under fairly general conditions. 

The following theorem, of the same type as (22.4), may 
often be used to determine the behaviour of a function for 
large values of the time : 

THEOREM VIII. // y(p) =f(p}jg(p\ where f(p) and 
g(p) are polynomials in p, the degree of f(p) being less than 
that of g(p), and if the roots of g(p) o are either zero or 
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have negative real parts * then if y(t] is the function whose 
transform is y(p) 

limy(t) = limpy(p) . . (22.5) 

t--co p*Q 

and, under the same conditions, 

poo 

y(t}dt = lim y(p) , . (22.6) 

J p-+Q 

The result follows immediately from the formulae of 
2 for expressing y(p) in partial fractions. An application 
of (22.6) has been given in (n.8). 

23. Network theorems 

Many of the general theorems of steady state network 
theory have immediate generalizations in terms of Laplace 
transforms. 

Here we mention only TheVenin's theorem because of 
its importance in simplifying the algebra in complicated 
networks. A special case of this theorem, which is ade- 
quate for many purposes, may be stated as follows : 

Suppose AB and CD are two pairs of terminals of a 
network ; let Z be the generalized impedance of the net- 
work looking in at CD with AB short-circuited ; also let 
v be the transform of the open -circuit voltage drop across 
CD due to voltage V() applied for t > o at AB with zero 
initial conditions. Then the transform of the current I 
in an impedance # connected across CD, due to voltage 
V(;f) applied at AB for t > o with zero initial conditions, 



s 



As an example, suppose the network is a filter circuit 
of in, T sections [Fig. 18 (a)]. In this case the quantity 
Z denned above is, by (16.10) with r = o, 

Z = % sinh 9 tanh mO . . (23.2) 

* This Theorern_ is frequently stated in the literature with 
no conditions on y(p~). Clearly it does not hold, for example, 
if y(t) = sin ait. 
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and the quantity v defined above is, by (16.11) with 
r = m i, 

v = V sech mO . . (23.3) 

Thus, by (23.1), the current I in an output impedance 
SQ is given by 



2 sinh 9 sinh m6 + # cosh tnff 

This agrees with (16-7), but the calculation is simpler, 
since The'venin's theorem allows the result for any terminal 
impedance to be deduced from the solutions for open and 
short-circuited terminations and these are often known 
and in any case are easy to calculate ab initio. 

24. The superposition * theorem 

THEOREM IX. If y^p} is the transform of y^t), and 
y t (p) is the transform of y z (t}, then yi(p)y z (p) the trans- 
form of 

ft n( 

(24.1) 

The proof f is not particularly difficult, but is omitted 
here as it depends on the transformation of double integrals, 
and involves only routine pure mathematics. 

* Also referred to as the Composition, Convolution, or Faltung 
Theorem, and sometimes as Borel's or as Duhamel's Theorem. 

t C. and J., 33 ; Doetsch, loc. cit., p. 161. The essential 
steps, which require justification and some conditions on the 
functions, are as follows : 



f-**yjfyi)du\ 
o Jo 



= f ^e- 



where the double integral is taken over the quadrant u > o, 
v > o. In this put u + v t, v = t', then t f can take all values 
up to t, and t all values up to ro, and the double integral becomes 

f "-' 
Jo 
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As a first application consider the differential equation 



to be solved with y == Dy = o , when t = o. The sub- 
sidiary equation gives 





Clearly, any system of ordinary linear differential 
equations with constant coefficients and with arbitrary 
functions on the right hand side can be treated in the same 
way. 

From the present point of view, the chief use or 
Theorem IX is the following : suppose we require the 
current I at some place in a network due to voltage 
applied at some point with zero initial conditions.* In 
usual way we find 



_ 

where z(p) is a known function of p. Then if x(i) is the 
function, found in the usual way, whose transform is 
!/#()> Theorem IX gives 

I = fv(*')*(* - t'}dt' = f V(t - t'}x(t'}dt' (24-3) 
Jo Jo 



Thus for any applied voltage V, which may be so com- 
plicated that we cannot evaluate its transform, or which 
may be given numerically, for example by an oscillogram, 
we can always express the current by a definite integral 
which may be evaluated numerically. In the same way, 
if the current in any part of a circuit is known, the current 
in any other part can be determined. 

* For non-zero initial conditions there will be additional terms 
on the right of (24.2) which can be dealt with as in Chapter I. 
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For example, for voltage V applied to an L, R, C circuit 
with zero initial conditions, we have as in (9.2) 

T 

and thus, if 2 > o, 

I = 4" f V(* - *'){ cos *' - /it sin *'}*-**' <fc' (24.4) 

W-LJo 

The function #() used above whose transform was 
i/#Q>) is, by (24.2), the current due to a unit impulsive 
applied voltage, V = 8(), V(p) = i. This is an interesting 
transient in its own right, and the present result shows 
that if it is known, the current due to any applied voltage, 
V, may be evaluated as a definite integral. The other 
simple and interesting transient in the circuit is the 
current Ij due to unit applied voltage, V = I, which, by 
(24.2) has the transform 



If IJL has been found as a function of t from (24.5), we 
may express the current due to any applied voltage V as 
an integral involving I 1} by writing (24.2) in the form 

' ' (34 ' 6) 



Now, by Theorem IX, the function whose transform is 

is 



&' (24-7) 

Then it follows from (24.6) and Theorem III that 

- t'}dt' y . . (24.8) 



since (24,7) vanishes when t o. Clearly various trans- 
formations of (24.8) are possible. 

Finally, it should be remarked that the importance of 
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Theorem IX is in the study of difficult or arbitrary func- 
tions, or functions given numerically by oscillograms ; if 
V is a simple exponential or trigonometric function, it is 
simpler to evaluate I directly by the methods used earlier 
than to evaluate the integrals (24.3) or (24.8). 

25. The inversion theorem 

This is a general theorem* which gives y(t) from y(p}> 
subject, of course, to conditions on y(t] or y(p}. Omitting 
these conditions, it states that if 



y(P) = f 

J 



o, . (25-1) 

o 

j ry+tco 

y(t)=-.\ eKy(X)dX, y > o. . (25-2) 

Z7TZJ y_oo 

In (25.2), A is a complex variable, 'and the integral is to 
be interpreted and evaluated by the theory of functions 
of a complex variable. If this theory, and in particular 
the calculus of residues, is known, (25.2) provides an 
important method for finding the function which has a 
given transform. For the problems of Chapters I and II 
which led to ordinary linear differential equations, this 
gives an alternative method for finding y(f), which, how- 
ever, is no shorter than those used earlier, so that for such 
problems there is no necessity to learn the technique of 
complex variable. For a thorough study of the applications 
of the method to partial differential equations, the use of 
the in version theorem is necessary. 

From the point of view of the historical development of 
the subject the correspondence between (25.1) and (25.2) 
is of great importance. Carson -f regarded (25.1) as an 
integral equation from which the solution of the problem 
could be obtained. Bromwich, and to a certain extent 

* C. and J., 29, 30. Doetsch, loc. cit., Kap. 6. The result 
is sometimes called the Fourier-Mellin Theorem. It is related 
to Fourier's and Mellin's inversion formulae, cf. Doetsch, loc. cit. 

f Electric Circuit Theory and Operational Calculus (McGraw 
Hill, 1926). 
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Wagner, Jeffreys, and McLachlan, regard an -integral of 
the type (25.2), with y(X) found from the usual subsidiary 
equation, as the solution of the problem. The Laplace 
transformation and its inversion theorem provides the 
link between these two ideas ; this point of view is largely 
due to Doetsch. 

26. The connection with Fourier's integral theorem and 

Parseval's theorem 

Fourier's integral theorem* states that, subject to con- 
ditions on y(t), if oj is real and 

TOO 

e imt y(t)dt, . . . (26.1) 



ro 

= 

j 



then y (t) = -L e- i(ot Y(a>)daj. . (26.2) 

277-J_co 

The simplest of the essential conditions on y(t) is that 

pCQ , 

\ 1X01^ be convergent. . (26.3) 

J-OO 

Parseval's theoremf states that, roughly provided these 
integrals exist, 

dt^ r |Y(o>)i 2 ^o>. . (26.4) 

ZTTJ _oo 

For comparison with the problems previously discussed, 
suppose that y(f) is real, and is zero J when t < o. Then 
(26.1) and (26.2) become : 

* For a proof see Carslaw, Fourier Series and Integrals, edn. 3, 
119 ; the result (26.2), (26.1) follows from his on using the 
exponential form for the cosine. See also Titchmar'sh, Theory 
of Fourier Integrals, i.i, 1.2. 

t Titchmarsh, loc. cit., 2.1. The theorem is often referred 
to as Rayleigh's theorem, and its use in physical applications is 
due to him [Phil. Mag., 27 (1889), 466]. 

J It is frequently stated that the Laplace transformation method 
contemplates only functions of this type. In the method as 
developed in Chapters I and II the problem was to find the solu- 
tion of a differential equation for t > o, which satisfied certain 
conditions at t = o, and in fact the solutions found would also hold 
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if 

then* -Rfi-**Ya)Xa>=XO> i 

(26.6) 

and iR^'YOXo. = o, if * 

f Jo 

provided | JXOI* is convergent. 

- Jo 

To -deduce (26.6) from (26.2) in this case, we notice 
that the right-hand side of the latter is equal to 



and, since y(t) is real, Y( CD) is the conjugate of Y(OJ), and 
the result (26.6) follows. 

Also, in this case (26.4) becomes 



Y(u) in (26.5) is called the Fourier transform of y(t). 
Since (26.5) is precisely the same as the equation defining 
the Laplace transform y(p) of y,(t), with p replaced by 
( z'o>), it follows that, if Y(w) exists, 

Y(w)=.y(-*co), . . (26.9) 

Thus the Fourier transform is the degenerate form of 
the Laplace transform if the latter exists when p is pure 
imaginary, also (26.6) is the degenerate form of (25.2) if 

for t< o; for example, the solution of (D + i)>" o with 
y i when t = o is found to be y = e~* for either t > o or 
t < o, though we are only interested in the former case. In 
this section and in 25 there is a difference in principle (though 
not in the results obtained for t > o) since a solution, la sought 
which is zero for t < o ; for example the solution of the differ- 
ential equation above, found from (25.2), would be y(t) ** o, 
t < o ; y(t) = i, t = o ; y(t) = e"*, t > o. 
* R is written for " the real part of ". 
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y can be taken zero. The Laplace transform exists for 
many functions for which the Fourier transform does not 
(for example, I, e* or sin f) but, if the latter does exist, it 
can be found from (26.9) and previous results. 

For example,* if y(t) = e~ at sm fit . . (26.10) 
so that 



_-_, 



. (26.Il) 

In this case (26.8) gives 

,. (26.12) 

^ v ' 



From the point of view of electric circuit theory the 
great importance of Fourier's integral theorem lies in its 
physical interpretation. The Laplace transform has no 
obvious physical significance ; it is a purely mathematical 
device for finding a solution this is not altogether a 
disadvantage as it implies that any difficulties can be form- 
ulated and treated mathematically, while it was on points 
of interpretation that the old operational method was 
obscure. On the other hand (26.6) expresses y(t] as a 
combination (that is the integral) of harmonic terms e~ imt , 
each with its own amplitude Y(o>)Ja), in the same way 
that a periodic quantity is expressed by Fourier's theorem 
as a sum of fundamental and harmonics, each with its 
own amplitude. 

Now suppose a voltage y(f) is applied to a circuit with 
zero initial conditions, and the current I in any branch of 
it is calculated, then, proceeding as in Chapter II we find 



* Fourier transforms of many functions of this type, with 
numerical illustrations and applications, are given in Burch and 
Bloemsma, Phil. Mag. (6), 49 (1935), 480. 
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where z(p) is known and depends on the circuit. 'Then 
if y(t] has a Fourier transform, this will be j7( z'o>), by 
(26.9), and (26.13) gives for the Fourier transform of l(t] 



Thus, by (26-6) 

-, ,>o. . (.6,5) 



This expression has a simple physical interpretation, 
since e- imt ]z( z'co) is the steady state current in the circuit 
due to voltage e~ i(at . Thus (26.15) may be interpreted 
by the statement that the current due to voltage y(t) is 
obtained by superposing the effects of all the component 
frequencies of y(t), each with its appropriate amplitude 
y( ia>)d<a. That is, we have an extension of the results 
of steady state alternating current theory to non-periodic 
applied voltages, although, for the solution of specific 
problems, this is inferior to the Laplace transformation 
methods described earlier. 

One application of (26.14), however, is new and of 
great importance. Applying Parseval's theorem (26.8) to 
(26.14) gives 



dot. . (26.16) 

TTJO "" <""' V ' 

The quantity on the left gives the energy absorbed in 
the branch of the circuit, alternatively it is the total re- 
sponse of a square law detector, which is of considerable 
practical interest, and (26-16) gives this as a real infinite 
integral without having actually to evaluate the current. 

For example, if voltage e~ a< is applied at t = o to an 
inductive resistance R, L we have 



o | (R Lzo>)(a ia>) | 2 
03 dw 

L 2 (t0 2 -f jB a )(w a + (X 2 ) [continued overleaf 
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_! _! 

~a tan ~ 'a Jo 



where j8 = R/L. 

In more complicated problems the integral (26.16) may 
be evaluated numerically, and approximations can be 
introduced, for example, with a selective circuit, \z(ioj}\ is 
very large except in a certain range of values of cu, and 
(26.16) need be evaluated only over this range. 



i. If y(f) is i, a, 3,.. .. in the regions o < t < T, T < t < aT, 
2T < t < 3T, and so on (a " staircase " function) 



2. If y(t) = i, o < t < tt ; y(t) o, t a < t < T ; and this 
is repeated with period T (a repeated pulse, or the difference of 
two staircase functions) 



3. If y(f) = | sin cat |, (full wave rectified alternating voltage) 

y(p} = ,^ ? ( 1,' 



4. A uniform beam is clamped horizontally at both ends. It 
carries a distributed load w (constant) per unit length in o < x < JZ 
and zero in $1 < x < /. Show that the deflection at any point is 



24EI 

5. If the square-wave voltage of 21, Ex. 3, is applied at 
f = o to an L, R, C circuit with zero initial current and charge, 
show that the current in the circuit at time t = asT + t' t where 
s is an integer and t' < T, is 
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using the notation (9.3) and assuming n z > o. Deduce that 
in the steady state due to a square wave voltage, the current at 
time t ' after the beginning of a cycle is 

2.E e-tf' sin nt' e-f(T+t') s j n n(T t') 
nL i + ze-pT cos nT + e~W 

6. In the problem of n, Ex. 2, show that for small values 
of the time 



I 



_ _ E / _ nW *^ \ 

~~ LV 6 + 60 ' ' 7 ' 



7. In the circuit of Fig. 7, an unknown voltage is applied at 
S at f = o with zero initial currents and charges. If the voltage 
drop across the first condenser is *;(), show that the voltage drop 
across the second condenser is 



r* 

n \ v(t') sin n(t t'}dt'. 
JQ 



8. If y(i) = (e~** e'P), t > o, show that its Fourier trans- 
form Y() is 



(a u)(/J z'o))' 

Show that the total energy absorbed by an inductive resistance 
L/, R from this voltage is 

RQ3 - q)'(q + jg + ^) 



a)' 
where p. R/L. 



CHAPTER IV 

PARTIAL DIFFERENTIAL EQUATIONS 

27. Introductory 

In the preceding chapters it has been found that problems 
leading to ordinary linear differential equations could be 
solved with only the elementary apparatus of 1-3- ^ n 
the study of partial differential equations it will appear that 
it is still easy to find the Laplace transform of the solution, 
but, since this usually involves more complicated functions 
of p than those considered in Chapter I, it is rather more 
difficult to derive the solution from it. The, Table of 
Transforms can of course be enlarged, but the most power- 
ful general method is the application of the inversion 
theorem, 25, combined with the use of contour integration 
to evaluate the complex integral which occurs in (z$-2)> 
For a complete study of the subject, particularly in its 
theoretical aspects, a knowledge of this technique is essential ; 
but it is possible to go far without it, for example most 
results on finite transmission lines, and some on infinite 
lines, can be obtained, in some cases rigorously, and in 
others with the proviso that for a complete treatment the 
complex variable would have to be used, although the 
actual process of writing down the solution would be 
precisely the same as that used here. 

The differential equations for a uniform electric trans- 
mission line will be found in 28, and, for shortness, only 
these will be studied, although many of the problems 
considered are problems in wave motion or conduction of 
heat in one dimension, with a slightly altered notation. 

28. The differential equations of the uniform transmission 
line 

We suppose the line to have resistance R, inductance L, 
capacity C, and leakage conductance G, per unit length. 
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Also let I be the current in the direction of increasing x 
at the point x of the line, and let V be the voltage 
drop across the line at this point. I and V will be func- 
tions of both distance x along the line and time t ; to 
emphasize this they will sometimes be written I(x, f) and 
V(x ft 

V ^"V, lj, 

To find the differential equations satisfied by the current 
and voltage in the line, consider the points A at x, and 
B at x -|- Sx. The current and voltage at A will be I and V, 

and at B they will be I + $x and V + -Sx, respectively. 
The portion of the line between A and B can be represented 

LSx RSx 
A H , > T 

A->. . T37T "W^ ""-^~ T "---- 

I 




X K *.fSoc 

FIG. 28. 

by the equivalent circuit of Fig. 28. Here the voltage 
drop from A to B over the series portion of the line is 

-- Bx, and the current in it is I + -Sx ; neglecting 
dtf d# 

terms in (S#) 2 this gives 



Also the voltage drop over the parallel portion HK of 

the line is V, and the current in it is S#, thus 

ox 



102 THE LAPLACE TRANSFORMATION 

From '(28.1) and (2,B\2) it follows that the differential 
equations for I and V are 






These are a pair of simultaneous partial differential 
equations for I and V. We wish to solve them for t > o 
with given initial values of I and V, which we shall write 
t(x) and V(x) to emphasize that they are functions of x. 
Explicitly these initial conditions are 

>- 1(#), as t -> o for fixed x, . (28.5) 
> V(#), as i ( -> o for fixed x. . (38.6) 

There are also boundary conditions to be satisfied at 
the ends of the line; these depend on the particular 
problem. For example, if the line is short-circuited at 
* = a, we must have V = o for x a and all t > o ; 
if there is open circuit at = a, the condition is I = o for 
x = a and all t > o ; if the voltage is prescribed at x = a, 
V(#, t) must be this prescribed function of t for t > o ; 
if the line is joined to some other impedance at x = a, 
there will be a condition involving both the current and 
voltage at x = a. 

29. The solution of partial differential equations by the 
Laplace transformation 

We consider first the system of equations (28.3), 
(28.4) which have to be solved for t > o with initial con- 
ditions (28.5) and (28. 6), 'and with boundary conditions 
to be specified. Proceeding * precisely as in Chapter I, 
4, we multiply the equations by e~^ and integrate with 
respect to t from o to oo. Then, just as in the proof of 
Theorem III, 

* As iii Chapter I assumptions about the function I are implied, 
and, strictly, the solution should be verified. 
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vWx t)dt 



u=o 
=*-!(*)+ pi . . . (39.1) 

where I is now a function of x and jp, and, when we wish 
to emphasize this, it is written I(x, p). 



where we have assumed that the orders of integration, a'nd 
of differentiation with respect to the parameter x, may be 
interchanged. 

Using these results, and the corresponding ones for V, 
we obtain from (28.3) and (28.4) and the initial conditions, 
the " subsidiary equations " 

(Lp + R)I(*. p) = - ^ + Ll(*) (29.3) 




(Cp + G)?(*. p) = - + CVW (29.4) 



Eliminating I gives the ordinary differential equation 
for?, 



where q* = (Lp + R)(C#> + G) . . (29.6) 

When V has been found, 1 is given by (29.3), that is, by 

1 = ~ Lp + R"^ + Lp+.R' ' 

It appears that the effect of the Laplace transformation 
with respect to t is to reduce a partial differential equation 
in x and t to an ordinary differential equation in x. This 
lias to be solved with boundary conditions which are the 
Laplace transforms of those originally prescribed they are 
"best written down as they arise in specific problems. 

The quantity q which appears in (29.6) is, in the general 
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case, the square root of a quadratic function of p, but in 
three important cases it takes a simpler form. These are : 
(i) The " lossless line ", R = G = o, for which 

q = p]c, where c = (LC)-* . . (29.8) 
(ii) Heaviside's " distortionless line " in which 
R_G_ 
L ~ C ~ p ' 

and thus - q = ( p -f- p)]c, where c = (LC)~*. . (29.9) 
(iii) The ideal submarine cable, L = G == o, for which 

q = -\/(RCp). . . (29.10) 

Another partial differential equation of importance is 
the equation of linear flow of heat (v temperature, K 
diffusivity) 



to be solved with v = v(x), when t o. The subsidiary 
equation, found as above, is 

d 2 v p_ i., , , N 

-j- 9 -v -- v(x) . . . (29.12) 

dx Z K K ^ ' ^ * ' 

and thus is of the same type as (29.5) in the case (29.10) 
of L = G = o. 

The equation of wave motion hi one dimension is 



where y is the displacement. This has to be solved with 
y = y( x \ ^ = Y(), when t = o. 

The subsidiary equation is 

d z v i) z -b T . 

2-=-)-), . (2 9 .14) 
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which is of the same type as (29.5) for the case (29.8) of 
R = G = 0. 

When the subsidiary equation (29.5) with its boundary 
conditions has been solved,* we know the transform 
*V(x, f) of the solution V(x, t). For problems on semi- 
infinite lines we rely on the Table of Transforms to find 
the solution from its transform. For problems on finite 
lines two methods are available : (i) an extension of the 
expansion theorem (2.8), which gives the solution in the 
form of a trigonometric series ; or (ii) an expansion in 
terms of the solutions for the semi-infinite line, which has 
a physical interpretation in terms of direct and multiply 
reflected waves. Both these methods will be used, treat- 
ing the latter first. 

30. The semi-infinite line 

We consider the problem of the semi-infinite line, 
x > o, with zero initial current and charge. First, suppose 
the end x = o to be maintained at constant voltage E for 
t > o. 

In this case the subsidiary equation (29.5) becomes 



___ = , - - 

The boundary conditions to be satisfied at the ends of the 
line are 

V = E, * = o, t > o, . . (30.2) 

V to be finite as x -> oo . . (30.3) 

The transforms of (30.2) and (30.3) are 

V=E/p, x = o, . . . (304) 
V to be finite as x -> oo. . . (30.5) 

* We could, as in 6, solve it by the ( Laplace transformation 
method, but, for the simple problems arising here, it is just as 
easy to write down the particular integral and complementary 
function in the usual way. 
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The general solution of (30.1) is 

V = Ae-<* + Be, 

and (30.5) requires B = o, and then (30.4) gives A = E/p. 
Thus, using the value (29.6) of q, 
p -p 



In the general case the transform involved in (30.6) is 
quite different from anything previously encountered, and 
in fact V turns out, after a difficult calculation which will 
not be given here, to be a complicated integral involving 
Bessel functions (cf. 35). In the important special cases 
mentioned in 29, however, V can be evaluated simply. 
For the lossless line, R = G = o (30.6) becomes 

p 
V = -e-vl6, t . . (30.7) 

P 
where c = (LC)-*. . . . (30.8) 

Then it follows from (21.1), with y(t) = E, that 

V:=EH[* -(*/<:)], . . (30.9) 

that is : V is zero up to the time x/c, at which a wave 
travelling with velocity c from the origin would arrive at 
the point x, and V has the constant value E at this point 
for greater times. 

For the distortionless line (29.9) we have 

p 
^ = e-(:P-t-P)/fl } . . . (30.10) 

p 

and' V = Eg-"/H(* xjc) . . '(30.11) 

In this case the disturbance is still propagated with 
velocity c, but there is attenuation as we move along the 
line. 

Some problems on the submarine cable, L = G == o, 
will be studied in 35. 

Suppose, now, that instead of applying a constant 
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voltage E at x o, we apply a sinusoidal voltage E sin cot. 
The only modification is that in place of (30.6) we have 



In the case of the distortionless line, q (p + p)/c, 
this becomes 



and V = Ee-p*>/ sin to(t - -W* - -) . (30.14) 

\ c / \ c / 

Still more generally, for any applied voltage, /(), we 
have in place of (30.6) 

V =/(>>-. . . (30.15) 

For the distortionless line this becomes 

17 =^)e-(i'+p)/c ) . (30.16) 

and, again using (21.1), 

. (30.17) 



This may be interpreted by the statement that the 
voltage at x is zero up to time x/c, and subsequently follows 
that at x = o with a time lag of x/c, and a reduction in 
magnitude by the factor e~i> a l c . 

In the next section the application to problems on finite 
lines of the solutions found above will be studied. 

31, The travelling wave solutions for the finite transmission 
line 

As an example of the method we consider the problem 
of the finite transmission line o < x <l, with zero initial 
current and charge. The end * = o is to be earthed, and 
the end x = / maintained at constant voltage E for t > o. 
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The subsidiary equation (29-5) becomes in this case 
^_ ? 2 V = o, o<x<L ' . (3 

The boundary conditions are 

V = E, * = /, t>o, . . ( 
V = o, * = o, t>o. . . ( 

The transforms of (31.2) and (31.3) are 



The solution of (31.1) which satisfies (31.4) and (31.5) is 

y = E m sinh g* . . (31.6) 

p ' sinh ql ' 

Here q is denned by (29.6) so that V is a complicated 
function of p. One method of proceeding is suggested 
by the device of expanding in a series of negative ex- 
ponentials which was used in 21, Example 6. We write 
(31.6) in a form involving negative exponentials, and 
expand the denominator by the binomial theorem. This 
gives 

E - 

J ' 

T? 

_ e -e(z-*)( I _ g-2a)(i -j- 



In the special case of the lossless line, for which, by 
(29.8), q = p/c, (31.7) becomes 

E r "\ ' 

"^ = J e -V(l-x)lc e -j>(l+as)/c _|_ e -D(3Z-aO/o _ ^ < I. ^1 .8) 
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And, using (21.1) 
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The graph of the function (31.9) is shown in Fig. 29. 
The voltage at x is zero up to time (/ x)/c, at which a 
wave travelling direct from the end x = I would reach the 
point x. The voltage then has the constant value E up 
to the time (I + x)/c, at which a wave travelling from the 
end x I and reflected back from the end x = o would 
arrive. From this time up to the time of arrival of a 
twice reflected wave, it has the value zero, and so on. 



Dir 
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(e-jcvfc ((.+x)/c (stbcyc (3Ux)/c 



FIG. 29. 

Returning to the general expression (31.7), it is seen 
that the terms of this are all of the form (i/p)e- 5 , which 
.appeared in (30.6) in the corresponding problem for the 
semi-infinite line. Thus, if the solution for this case is 
known, the solution for a finite line can be expressed in 
terms of it. 

As a second example, we consider a lossless line of 
length /, initially charged to constant voltage E. At 
t = o the end x = o is earthed, the end x = I being left 
insulated. 

In this case the subsidiary equation (29.5) becomes, 
writing c = (LC)~*, 
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which has to be solved with 

?=o, a = o, . . (31-") 
I = o, x = L . . (31.12) 

Also it follows from (29.7) that 

1 'f. . . ( 3 ,r 3 ) 

Lpdx 

The solution of (31.10) which satisfies (31.11) and (3 1. 1 2) is 

, v 

p cos ' ' ' ' u 

_Ef g-J"/c[ I 

= j>V i. 

p 
= {l - e-pa>/C 

Thu, ^j-Hfi-^- 

E \ c) c 



Also by (31.13) 



k 



.:. (31.15) 






The most interesting quantity is the current at = o. 
At this point 

E 



(3I- 1 ?) 
that is, a square wave of amplitude E-\/(C/L) and period 
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32. The finite line with terminal impedances 

Suppose that at x = I, where the voltage is V z and the 
current is I r , the line is terminated by impedance z^ 
Then, for zero initial conditions in this impedance, the 
boundary condition at x = I ^..Q y.* 

is just the subsidiary equation y ... 

for the impedance z^ namely Al^ 




Similarly, if voltage v is 
applied through impedance 
Z Q to the line at x = o ;, 
writing V for the voltage at IGl 3 " 

the point x = o of the line and I for the current there, 
the boundary condition at * == o, for zero initial con- 
ditions in the impedance % Q is 

As an example, suppose that a lossless line of length / 
is initially charged to voltage E, and that at t o it is 
discharged through a resistance R at x = o, the end 
x = I being insulated. In this case the*subsidiary equation 
(29.5) becomes 

where c = (LC)~*. At x /, we have 

I = o, x = 1, . . . (32.4) 
and at x = o the boundary condition (32.2) is 

R I = -?.. . . (32.5) 
Also we know from (29.7) that 

1=-^^, .. .- (32-6) 
and thus (32.5) becomes 



= 0. 
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The general solution of (32.3) which satisfies (32.4) is 

V = ? + A cosh p(l-x)lc. . (32.8) 

Substituting this in (32:7) gives 

.R . . pi A , pi , E 
A=4 sinh - + A cosh-i- H- T = o. 
Lc c c p 

Using this value of A in (32.8), we find 



- 
p R sinh (p//c) + Lc cosh (p// 

The current is obtained from (32.6) which gives 

/ 
' U ' 



jp{R sinh (pile) + R x cosl 

where R x = Lc = V(L/C). For shortness we evaluate 
only the current at x = o. In this case, expressing (3 2 ^) 
in negative exponentials, and writing a=(Ri I 
we find 

E(r -* 

. + RO)(I 

E 



(a 



Therefore 

E 



. 



so that the current has a series of constant values, diminish- 
ing in magnitude. If R = R x = -\/(L/C), (the character- 
istic impedance), there is a single pulse of current E/aR 
for time zljc, and zero current subsequently. 
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Clearly any terminal impedances can be discussed in 
the same way, but since the coefficients of the later ex- 
ponentials are usually rather complicated functions of p, 
the method is most suitable for use when the solution is 
needed for fairly small values of the time. 

33. The expansion formula for finite lines 

In chapters I and II the transforms of the solutions of 
problems on ordinary differential equations were usually 
found to be of the form 

<p(p) = ... r^^.i) 

where f(p) and g(p] were polynomials in p. The solution 
ff>(f) was found from (j>(p] by expressing (33.1) in partial 
fractions, or by the equivalent technique of Theorem II, 
for example, by the result (2.8) that if a^ . , ., n are the 
zeros ofg(p), and these are all different, then 



In the problems on finite lines in this chapter, rather 
more complicated transforms have appeared. The new 
feature of these is the occurrence of quotients of hyper- 
bolic functions of q. For example, the expression 

' (33-3) 



pcohpl/c 

appeared in (31.14). Since cosh % and sinh % can be 
represented by the infinite products * 



> (33-4) 

;nb) > (33-5) 



* Carslaw, Plane Trigonometry, edn. 3, 168. 

8 



a quotient of hyperbolic functions such as (33.3) may still 
be regarded as being of the type 

C33-) 

except that nov?f(p) and g(p) may have an infinite number 
of factors, instead of a finite number as in (33.1). Then 
if a lt a 2 , . . . are the zeros of g(p), and if these are all 
different, it is plausible to assume that (33.2) can be ex- 
tended to this case : that is, if $(p) is given by (33.6), 

- ' ' (33 ' 7) 



Just as in 2 (33.7) applies only to the case in which 
g(p) has no repeated zeros ; if it has, the corresponding 
extensions (2.19) or (2.20) must be used. 

The result (33.7) will be referred to as the expansion 
formula ; it has not been proved, and is merely suggested 
as a possible extension to the case of an infinite number of 
zeros, of the result (2.8) which was proved for the case 
of a finite number of zeros. - That is precisely the approach 
which was, used by Heaviside. The justification for the 
use of (33.7) is as follows : the rigorous pure mathematical 
method of finding </)(t) from (33.6) consists of the applica- 
tion of the inversion theorem (25.2), followed by contour 
integration and a certain amount of pure mathematical 
discussion, and, in the case of finite lines,* the final step in 
this procedure is equivalent to the use of (33.7). Thus it 
may be taken that there is a solid theoretical foundation 
for results on finite lines (or on wave motion, or conduction 
of heat, in a finite region) obtained by the use of (33.7) J 
this remark does not apply to semi-infinite lines,^ for example 
(33.7) must not be applied to (30.6). 

* Any applied voltages must be combinations of the functions 
of Table I, and the lines must not be terminated by their char- 
acteristic impedance. 

t Or to lines terminated by their characteristic impedance 
in .this case the transforms encountered are of the same type 
as those for semi-infinite lines. 
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As a first example of the use of (33.7), we apply it to 
find (j)(t) from (33.3). First the zeros of the denominator 
must be found. This could be done by putting * z = plfc 
in (33.4), but, from the point of view of the subsequent 
calculation, the following alternative method is better : 
clearly = o is a zero, and the others are 

(in + i We . 

=^ - ___ } n==0> I( 2,..,, (33:8) 

since, for these values of p, 

, pi . c(zn 4- i We l\ (zn 4- iW 

cosh = costu - - -~ . - V = cos v - ! = o. 

C \, 21 CJ 2 

To apply (33.7) to (33.3) we need 

d f* ^P l \ \.P l , 1 P \P l / \ 

-j- (p cosh )= cosh 4- smh . (33.9) 
dpY cJ c ' c c VJJ ^' 

Thus __(*cosh H = i. . (33.10) 

UA cjj f=0 ** ' 

Also, since cosh (pl/c] = o when p has one of the values 
(33.8), we find from (33.9) 

f d ( pl\~\ (zn rf i>'77 . (zn + i)zV 

-5- (p cosh ) =- - smh- - 

\_dp\ ^/Jj)=(2n+l)Mro/2Z 2 2 



Using these results, (33.7) gives from (33.3) 

Conjugate V 



I r COS 



L 2 

jc (2M+ I)TTC< . , 

- cos - - L y , (33. iz) 

aZ 

* (33.4) and (33.5) can be used to show the actual factors of 
numerator and denominator : this is often useful in cases where 
a factor can be cancelled between them, or where some exceptional 
factor arises. 
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where the term i has arisen from the zero p = o of the 
denominator, 
Using this result in (31.14) gives 

00 

. . 

(33 - I3) 



A solution of the same problem, expressed in terms of 
successively reflected waves, was given in (31.15). The 
connection between this and (33.13) is of course that at 
any point x, the solution (31.15) is a periodic function of 
t' with period 4//c, and if this function is expressed as a 
Fourier series, the result is (33.13). 

As another example of the application of (33.7) we deter- 
mine V for the first problem of 3 1 from its transform 
(31.6), namely 

, , 

' (33 ' I4) 



considering various forms of g-, which was defined in 
(29.6) by 

? 2 = (Lp + R)(Q> + G). . (33-15) 

The first point to notice is, that if we use the infinite 
product (33.5) for sinh % in (33.14), we get 

E gl + Vl 'g V . 



. . , 

. . t33 ^ 

It follows that although V, as given in (33.14), appears 
to be a function of q (which is the square root of a function 
of p), it in fact involves only q 2 , and thus is a function of p 
which does not involve square roots. The same result 
holds for all transforms which arise in problems on finite 
lines. 

The general case of (33.14) will be studied in Example 3 
below ; special cases are given in Examples i and 2. 
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EXAMPLE i. The '" lossless " line, g = pjc, where 
c = (LC)-i 
In this case (33.14) becomes 

. E sinh pxlc . . 

V = . , ., -. . . (33.17) 
p sinh pi c 

One zero of the denominator is p = o, and the others are 
p=inirc[l, ni, 2... . (33.18) 

. , fimrc 1} . . 

since sinh < = . - V = t sm mr = o. 

{I c) 



AI T. ^ 

Also -j- p sinh == - cosh - 

|_0p\ C/Jp=iTaa>ll LC Cjj, iwtc // 

= (} n imT, n = i, 2, ... (33.19) 

The zero j> = o of the denominator is clearly exceptional 
since the numerator of (33.17) also vanishes for this value. 
The simplest way to see what happens in such a case, and 
to evaluate the corresponding term in the expansion formula, 
is to use the series for sinh z and cosh z, namely 

% 3 z 5 

= s +-J+TJ + -. (33-20) 



(33- 21 ) 
ctor. Using 

+ (pWI^} + . . .] 



in (33.17), and cancel any common factor. Using (33.20), 
(33.17) becomes 






and, in this form, the numerator of V does not vanish when 
p = o, and the denominator has only a single zero there. 



n8 



Using this result for p = o, and (33.19) for the zeros 
(33.18), we get finally from (33.17) and the expansion 
formula (33.7) UJ " F 



F* 

V .-: 



-- __ sn . e nxo __ Conjugatej 






. 

(33-24) 



The minor difficulty which was encountered above at 
p = o arises in many problems, and it is always advisable 
to see precisely what happens at such a point by substituting 
the series (33.20), (33.21), or the products (33.4), (33.5) for 
the hyperbolic functions, and cancelling any common 
tactor. In calculating the contribution from such a point, 
instead of using a term of the type [f(a r ]e a ^/g'(a r y\ which 
occurred in (33.7) and (2.8), one of the equivalent type 
(2.10), namely 



could have been used. In the problem above, this gives 



which is the first term of (33.24). 

Note also that if, after cancelling any common factors 
as in (33.22), the denominator of V has a multiple zero, 
a term of the type (2.20) must be used for this zero. 

EXAMPLE 2. The submarine cable, q = (p/ic)*, where 
K i /RC. 

In this case (33.14) becomes 



p*mhlV(PlK)' ' ' '' 

As in Example i, p = o is an exceptional value at which 
both numerator and denominator vanish. Usin .20 



we find 
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and 



The other zeros of the denominator of (33.26) are 

p = - KM// a , n = i, 2, 3 . . ., (33-28) 



since sinh l( ^- ) = db sinh imr o. 

Also \~ 



. (33.39) 

-KnWIP 

Thus we have from (33.26) and (33.7), using (33.27) and 
(33.39) 



v = +f5y a *w| 

(33-30) 



EXAMPLE 3. The general case of (33.14) with q given by 
(33.15). Here p = o is a zero of the denominator, and 

r^-(psinh?oT = [sinh gZ]j, = o = sinh Z(RG)i (33-30 

\_dp Jcp=Q 

The other zeros of the denominator are found from the 
values of q which make sinh ql zero ; these are q = o and 

q = -j- //, = i, a . . ., (33-3 2 ) 
The zero 9 = need not be considered, since a factor 
q has cancelled between numerator and denominator in 
(33.16). Using the value (33.15) of q, (33.32) gives 



(Lp + R)(Cp + G) + f- = o, n = i, 2 . . ., (33-33) 
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This is a quadratic which gives two values of p correspond- 
ing to each value of n ; it could also have been deduced 
from (33.16). Using the notation 

P = (R/aL) + (G/aC), a = (R/aL) - (G/aC), 

c = (LC)-*, . . (33-34) 
(33-33) becomes 

p 2 + 2pp + p 2 o- 2 + ^ =o, n 1,2 (33-35) 
The roots of this are 

-Piv n , . (33-,3 6 ) 
where v n = [(wW/Z 2 ) - a 2 ]*, . . (33-3?) 

and we suppose for simplicity that v n is feal for all values 
of n, that is (TTC[I) Z > a 2 ; if this is not the case (33.35) 
will have real negative roots for some low values of n, and 
these must be treated separately. 
Si As usual we need 



-j-(p sinh ql = \lp- cosh ql 

Lfy *J P =- P +. n L'dp tf J, P4 * 



cos 



W!s . (33.38) 



Using (33.31), we have finally by (33.7) 

_ T Esinh*(RG)* , -n-Ec 2 . 
V = _ x ' . -\ -- p-pt 
sinh/(RG)i ^ Z 2 







Conj.j 



sinh 



" 



where ' tan S n = p/v n . . . (33-39) 
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The examples above seem sufficient indication of the 
method of applying the expansion formula to problems on 
finite lines. From the results (33.13), (33.24), (33.30), 
(33.39) it appears that the expansion formula gives the 
solution as an infinite series, whose terms are products of 
trigonometric terms in x and exponential or trigonometric 
terms in t. These are quite different in form to those 
obtained in 31, 32 in terms of multiply reflected waves. 
Each type of solution has its own advantages, so that it is 
hard to give any general indication of which method it is 
better to use preferably both should be tried in any 
specific problem. However, the position as seen from the 
results of the examples worked out above is as follows : 

(i) For the lossless and distcfrtionless lines with simple 
terminations, both methods are available, but the solution 
in terms of multiply reflected waves gives the actual wave 
form of the disturbance with a simple physical interpreta- 
ion, and thus is preferable. 

(ii) For more general lines the solution by the expansion 
formula is very little more complicated than that for the 
lossless line, while the solutions in terms, of multiply 
reflected waves involve much more complicated functions 
(not given here). 

.(iii) For lines with terminal impedances, the solution in 
terms of multiply reflected waves is useful in general only 
for the first few waves, that is, for not too large values of 
the time, while the solution obtained by the expansion 
formula (discussed in 34) is available for all values of the 
time, though it is often difficult to determine the zeros of 
the denominator of its transform. 

34. The application of the expansion formula to the finite 

line zoith terminal impedances. 

Problems of this type are treated precisely as in 33, 
the only new feature being that the transforms have more 
complicated denominators whose zeros may be difficult to 
find. As a specific example we consider a lossless line of 
length /, insulated at x = I, and suppose that at t = o a 
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condenser of capacity C , charged to voltage E, is discharged 
into the line at x o. 
The subsidiary equation for the line is, writing c = (LC)-*, 

<*</, . (34-1) 



= = O, X = 7 /I A -9^ 

dx ' ' ' ' \34 >2 -J 

If V and I are the voltage and current at the point 
* o ot the line, the subsidiary equation' for the condenser 
C is ^ 

_.JLi_ ? _E ( , 

c oP p' ' ^34'S; 

that is, using (29.7), * 

i dV ^L E 
T <n +.% -T- V , when x = o. . (34.4) 

-L/V^np Q.X i) \<J~ TV 

This is the boundary condition to be satisfied at * = o. 
ihe solution of (34.1) which satisfies (34.2) is 

V = A cosh [p(l - x)/c)]. 
Substituting this in (34.4) gives 

A (L(f 
Thus ^ 

V = 1^2 cosh j>(/ x)/c 

Cc ' sinn (pi/c) + (pC /Cc) cosh (pljc)' 
The zeros of the denominator of (34.5) arep = o, and 

P= ^oc n //, n= i, z . . ., . (34.6) 
where the n are the positive roots of 

sin a + ka. cos a = o, . . (34.7) 

where p /in / n \ 

*-C //C . . . (34.8) 

is the ratio of the terminal capacity C to the capacity of 



the whole line. Approximate values of these roots can 
be found from the intersections of the curves 

y = tan a j 
and y ka.f' 

as shown in Fig. 31, and these values can be improved by 
interpolation. 




,y^ 



** 



FIG. 31. 



To find V from (34.5) by the expansion formula (33.7), 
we need 



= (l/c)[k + i Hh A: 2 a] cos a, 
where k is given by (34.8). Also 



(34.9) 



. (34.-) 
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Using (34.9) and (34.10), the expansion formula gives 






cos 



35. Further discussion of semi-infinite lines 

In 30 the transforms of the solutions of problems on 
semi-infinite lines were found to contain the factor 

- . . . (35- J ) 
where q = */{(Lp + R)(Cj> + G)}, . (35- 2 ) 

and in 31, 32 the travelling wave solutions for finite 
lines were found to involve similar factors. It is thus 
very important to know the functions which correspond 
to transforms of these types. In the case of the lossless 
and distortionless lines these can be found by the use of 
Theorem VI ; in all other cases, either the inversion 
theorem, or a table of transforms, must be used. Perhaps 
the most useful collection of transforms is given in Campbell 
and Foster x * " Fourier Integrals for Practical Applications ' , 
in the relevant part of which, Part 8, most of the results f 
can be read as Laplace transforms, p being our p, and g 
being written for our t. Many results for the general 
case (35.2) are given in their Nos. 860.5 onwards, and, 
for the submarine cable case, L = G = o, earlier. 

Here we consider only the submarine cable case, which 
corresponds also to linear flow of heat, and give an ex- 
tension of our table of transforms which will cover a 
number of problems of this type. The simplest and most 
important results follow from the known definite integral 

f V&*u3-<a*/4w) fa = ?Le- a i>, . (35-3) 

JO 20 

*Bell System Technical Monograph, 8-584. See also 
McLachlan and Humbert, Formulaire pour le calcul symbolique 
(Gauthier-Villars, 1941). 

f Namely, those stated for g > o. 

j Gibson, Treatise on the Calculus, edn. 2, p. 470. 



PARTIAL DIFFERENTIAL EQUATIONS 125 

where a and b are positive. , Putting b z = p, and sub- 
stituting M = \/t in this, gives 

*-W^ = j-'* . (354) 

Integrating (35.4) with respect to the parameter * a, 
from a to oo, gives 



f 00 

I'-" 


t tdt^e *l*da p^e ** dt 


1 (35-5) 


Writing 


Qffx ~~ 1 c~~^ du 


(35-6) 


erfc 


2 f 00 

x = i erf* = vl e~ uz du, 

*)x 


(35-7) 



for the error function (numerical values of erf x are given 
in most books of tables) (35.5) becomes 

e -ai>* 

- (35 - 8) 



In Table 2 the results (354) and (35.8) are given, together 
with some others which are useful in more complicated 
problems. Other transforms can be deduced from these 
by Theorems' I, III, IV, IX, also by differentiating or 
integrating with respect to a parameter, just as (35.8) was 
deduced from (354). For example (35.9) follows from 
(35.10) by differentiating with respect to the parameter a 
(or alternatively it may be derived independently from 
(3S-3) by putting u t~*). Also (35.11) may be deduced 
from (35.9) by Theorem IV. The results (35.12)^0 
(35.14) may be proved by a rather more elaborate discussion 
of the same type. 

* The processes of differentiating or integrating under the 
integral sign with respect to a parameter are most useful, since 
they often allow a chain of transforms (such as (35.9) to (35.11), 
or (35.12) to (35-14)) to be deduced from a single transform. 
It must always be understood that the use of these processes 
needs justification, but this is usually readily supplied for integrals 
of the types considered here. 
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Finally, the Laplace transform of t v , where v is not 
restricted to be a positive integer as in Table I, must be 
mentioned. This is 



P Jo 

TABLE 2 
In this Table a and b are positive or zero 



y(P). 



y(t). 



(35-9) 



erfc (Jflfi) 



-^ + 



e al>+W er f c _)_ 



P(b 



. (35.13) 
(35.I4) 
. (35-15) 



If v> i, the integral in (35.16) exists and is denoted 
by P(v + i) ; this gives the result (35.15). Numerical 
values of the F function are given in most books of tables. 
Its most important properties are 

and r(ri) = (n i) !, if n is a positive integer, (35.17) 
so that (1.5) is a special case of (35.15), as it should be. 
Also 

so that the special case a = o of (35.10) agrees with (35.15)- 
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As examples of the use of Table 2 we consider the follow- 
ing problems on the semi-infinite submarine cable. 

EXAMPLE i. Constant voltage E applied at x = o at 
t = o with zero initial conditions. Here, from (30.6) with 
L = G = o, 

V=je-*(Pl"* . . (35.19) 

where K = i/RC. . . . (35.20) 

Then it follows from (35.11) that 



Using (29.7), the current I is to be found from 

T , i dV EC* . , .1 , . 

T = = g-^WK)* fac.22) 

R dx (Rp)* v 

Thus, from (35.10), 



If we had only required the current I at x =^o, (35.22) 
becomes in this case 

(35 ' 24) 



and, from (35.15) and (35.18), 

' ' ' (35 ' 2S) 



EXAMPLE 2. The problem of Example i, except that the 
voltage is applied to the line through an impedance * - 
As in 30, the subsidiary equation for the line is 

fj_V = o, *>o, '. . (35.26) 
ax* K 

where K is denned in (35.20). The solution of (35.26) 
which remains finite as x -> oo is 

V = Ae-V(*M. . . (35.27) 
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Also, by (29.7), 

The constant A is found from the boundary condition 

T E rr' 
** l =p-~> x = - 

That is, P. 

Thus V = ^ JVK " e-V(p/) . (35.29) 

If # is a resistance, R,,, it follows from (35.14) that 
V = E erfc -^ . Ee**+* a|rf erfcf ^ r -f k-\/(i<f) \ 

(35-30) 



where k = 

Again, if # is a condenser of capacity C , (35.13) gives 



v = 

(35-3 1 ) 

For other values of # , the expression i/feop^ + R*) /n 
(35.29) must be expressed in partial fractions whose 
denominators are linear in^>*. 

EXAMPLES ON CHAPTER IV 

1. Voltage E sin (<at + a) is applied at t = o through a re- 
sistance RQ to a semi-infinite lossless line x > o, with zero initial 
current and charge. Show that the voltage in the line is 

T7"D r*/ **\ i / **\ 

R ^7 R Bfc[ (* - *j +ajH (i - * J , 

where c = (LC)~*, Rj = V(L/C). 

2. Constant voltage E is applied at t = o at the end x = o 
of a lossless line of length /, the end x / being insulated. If 
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the initial current and charge are zero, show that the voltage 
at any point of the line is 

~ (zn + z)l - 



3. A lossless transmission line extends from x I to x = a, 
the end x = a being insulated, and the end / being earthed. 
At t o there is no current in the line, the portion o < x < a 
is charged to voltage E, and the portion I < x < o is uncharged. 
Show that the current at x == I at time t is 



zn(l + a} + 



-] 



(2tt + i)(Z + a) + a 

L G 

where c (LC)~*. 

4. Constant voltage E is applied at t o through an inductive 
resistance L , R to the end x o o'f a lossless line of length Z, 
the end x = I being insulated, and there being zero initial current 
and charge. Show that the current at w = o .is 



aERi 



where c = (LC)-i, R t = (L/C)*, and a = (R + RJLp. 

5. A submarine cable (L = G = o, and K= i/RC in (29.10)) 
of length I has zero initial current and charge. The end x = o 
is insulated, and a constant voltage E is applied at x = I, show 
that the voltage at any point is 



, 

7T 



6. A lossless line of length Z has zero initial current and charge, 
and is insulated at x = /. Voltage E sin wt is applied at x o 
for t > o, show that the voltage at any point is 

E cos [tu(Z x)jc] sin <of 



^- 

8E<0/C 2. 



_ (zn 
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' 3nd w is not to be e 1 ual to ( 2 



tor any . 

7- A submarine cable of length Z has the end x = I earthed, 
and constant voltage E is applied at x = o with zero initial con- 
ditions ._ bhow, using the method of 31, that the current at 
x o is 
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